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1.0  INTRODUCTION 

Neutron  riif fraction  was  first  demonstrated  in  1936  using  a  paraffin 
moderated  RaBe  neutron  source  and  magnesium  oxide  crystals.  Because  of 
the  lov;  source  strength,  no  quantitative  information  was  obtained.  Neu- 
tron diffraction  had  to  await  the  development  of  nuclear  reactors  and 
their  high  intensity  thermal  neutron  beams  before  becoming  a  valuable 
research  tool. 

The  first  neutron  spectrometer  (1,  2)  was  constructed  at  Argonne 
National  Laboratory  in  1945.  The  range  of  this  early  machine  was  0.04 
eV  to  65  eV.  Since  that  time,  the  design  of  neutron  spectrometers  has 
undergone  many  changes.  Many  of  the  modern  instruments  have  three  axes. 
They  are  thus  able  to  bombard  a  sample  with  a  monochromatic  neutron  beam 
and  then  analyze  the  diffraction  pattern  obtained  for  any  given  energy. 
Some  (3)  of  these  instruments  are  also  able  to  apply  a  magnetic  field 
to  the  sample  and  to  vary  its  temperature  from  near  absolute  zero  to  over 
2000°C.  Naturally,  data  is  obtained  completely  automatically. 

In  the  ensuing  years  since  1945,  the  lower  energy  limit  (4)  of  the 
instrument  has  been   decreased  to  around  .00062  eV.  This  has  been  done 
through  the  use  of  special  neutron  filters  to  reduce  higher  order  contam- 
ination of  the  incident  beam  and  by  the  use  of  special  crystals  and  liquid 
hydrogen  moderating  chambers. 

Little  success  has  been  achieved  in  extending  the  upper  limit,  how- 
ever. As  the  upper  limit  is  approached,  the  resolution  of  the  instrument 
greatly  decreases  and  the  diffracted  beam  intensity  becomes  very  weak. 


Generally,  a  neutron  chopper  is  preferred  for  measurements  of  neutrons 
having  energies  in  excess  of  10  eV. 

The  neutron  spectrometer  has  several  advantages  over  its  main  com- 
petitor, the  neutron  chopper.  The  neutron  spectrometer  costs  consider- 
aoly  less,  is  much  easier  to  operate,  requires  less  floor  space,  and  is 
much  less  dangerous  in  operation.  The  chopper,  however,  covers  a  much 
wider  energy  range  although  drastic  changes  in  energy  usually  necessitate 
changing  rotors.  Also,  the  duty  cycle  of  a  chopper  tends  to  be  very  low, 
thereby  reaucing  the  average  intensity  of  the  chopped  neutron  beam. 

In  the  past  several  years,  some  of  the  best  features  of  the  chopper 
and  the  neutron  spectrometer  have  been  combined.  A  small  chopper  (5)  is 
often  introduced  into  the  neutron  beam  prior  to  the  diffraction  crystal. 
The  resolution  of  the  chopper  is  purposely  made  very   poor  in  order  to 
provide  a  high  duty  cycle.  Its  function,  then,  is  to  provide  a  means  of 
removing  the  faster  neutrons  from  the  diffracted  beam  in  order  to  reduce 
its  harmonic  content.  Sometimes,  even  two  phased  choppers  are  used  for 
this  purpose.  When  operateo  for  this  purpose,  a  chopper  is  usually  re- 
ferred to  as  a  mechanical  monochromator. 

For  certain  types  (7)  of  experiments,  neutron  diffraction  methods  are 
preferable  to  those  of  x-ray  diffraction.  The  scattering  cross  sections 
for  x-rays  are  proportional  to  the  atomic  number  of  the  atom  doing  the 
scattering.  Therefore,  this  cross  section  is  low  for  the  light  elements 
such  as  hydrogen  and  carbon.  Also,  there  is  little  variation  between  two 
atoms  having  nearly  the  same  atomic  number.  Since  neutrons  are  scattered 
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It  has  been  noted  (6)  tnat  the  energy  stored  in  the  spinning  rotor  of  a 
fast  chopper  may  be  as  high  as  10^  foot-pounds. 


primarily  from  the  nucleus,  the  neutron  scattering  cross  section  may 
vary  considerably  between  different  isotopes  of  the  same  element.  Thus, 
neutron  diffraction  is  an  important  tool  in  structure  studies  of  materials 
containing  either  the  very   light  elements,  or  of  materials  containing 
elements  with  very   close  atomic  numbers. 

The  other  type  of  experiment  is  the  study  of  magnetic  materials. 
Since  the  neutron  has  a  magnetic  moment,  it  is  capable  of  providing  infor- 
mation about  the  atomic  changes  that  occur  when  materials  are  magnetized. 

Of  course  the  most  obvious  use  for  a  neutron  spectrometer  is  to  deter- 
mine the  neutron  spectrum  emerging  from  the  various  experimental  facilities 
at  a  nuclear  reactor.  One  may  then  study  the  effect  of  changing  the  reac- 
tor core  configurations  upon  the  emergent  neutron  spectrum.  This  informa- 
tion may  then  be  used  to  allow  the  neutron  spectrum  to  be  somewhat  tailored 
to  suit  other  types  of  experiments. 

Another  use  of  the  neutron  spectrometer  is  in  the  determination  of 
neutron  cross  sections  as  a  function  of  energy.  For  this  type  of  measure- 
ment, a  correlation  is  made  between  a  spectrum  taken  with  the  sample  in 
the  diffracted  beam  and  one  taken  without  the  sample. 

The  neutron  spectrometer  described  later  in  this  work  was  constructed 
with  the  possibility  that  it  could  be  used  as  a  piece  of  experimental  appa- 
ratus in  an  undergraduate  laboratory  course  at  the  senior  level.  The  fol- 
lowing theory  has  been  written  with  this  application  in  mind.  In  addition, 
complete  instructions  have  been  included  to  allow  proper  care  and  alignment 
of  the  instrument. 


2.0  THEORY 

2.1  Basic  Concepts 

A  neutron  spectrometer  is  the  neutron  analogue  of  the  familiar 
x-ray  spectrometer.  Each  instrument  has  a  source  of  radiation,  a  dif- 
fracting crystal,  and  a  detector.  The  main  difference  in  the  outward 
appearance  of  the  two  types  of  spectrometers  is  one  of  size.  The 
common  x-ray  spectrometer  may  usually  be  placed  upon  a  table  three 
feet  square  and  weighs  less  than  500  pounds  (not  counting  the  necessary 
electronic  equipment).  A  typical  neutron  spectrometer  may  utilize  a 
ten  foot  square  and  weigh  many  thousands  of  pounds.  The  important 
elements  of  a  neutron  spectrometer  are  shown  in  Figure  1. 
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Fig.  1.  Schematic  diagram  of  a  neutron  spectrometer 


Every  portion  of  the  system  indicated  in  Fig.  1  which  either  comes 
in  contact  with  or  generates  the  neutron  beam  that  emerges  from  the 


collimator  will  have  an  effect  upon  all  measurements  taken  with  the 
system.  In  order  to  better  understand  the  entire  process,  let  us 
first  look  at  the  phenomena  of  neutron  diffraction. 

For  nearly  four  decades  it  has  been   known  that  under  appropriate 
circumstances  behavior  of  the  elementary  nuclear  particles  could  be 
predicted  by   assuming  that  they  were  associated  with  a  wave  having 
a  wavelength 

where  x  is  the  particle  wavelength, 

h  is  Plank's  constant, 

m  is  the  particle  mass,  and 

v  is  the  particle  velocity.  v 

From  this  point  on,  it  will  be  assumed  that  equation  1  is  a  valid  descrip- 
tion of  the  neutron  during  the  process  of  neutron  diffraction. 

Fig.  2  shows  the  relationship  between  the  neutron  velocity,  wave- 
length, and  kinetic  energy,  E9  where  the  latter  is  given  by   equation  2, 

E  =  1/2  mv2.  (2) 

It  is  instructive  to  observe  the  ranges  of  the  various  quantities  plotted 
in  Fig.  2.  Neutrons  having  the  energies  indicated  are  found  in  great 
abundance  within  the  core  of  thermal  nuclear  reactors.  Furthermore,  the 
corresponding  wavelengths  are  within  the  order  of  magnitude  of  interplanar 
spacings  of  common  crystals.  Hence,  suitable  neutron  sources  and  diffracting 
materials  are  available. 


c  +->  E 
o  en  o 

S_    C  S- 

4->      O    +-» 

Z5  i —  l/> 

<1J     CJ  Ol 

Z   >  c 

fO  TO 


> 

S- 

o 
a 


00 

> 

> 

r- 

CJ 

-t-> 

C! 

c 

>> 

a 

o 

r— 

s_ 

a 

a 

> 

c 

a 

—■ 

■a 

c 

fD 

>> 


u 

c 

a 

> 

c 
c 
s- 

a 


csj 


o 


C 

o 

S- 

+-> 
a 


o 
o 


o  +->  u 
s-  t-  a 

•P     UU1 

3    O  "^ 
CD  r—    EI 


2.2  Derivation  of  Bragg  Law 

Consider  a  perfect  crystal  (8)  placed  in  a  parallel  beam  of 
monochromatic  neutrons.  The  angle  between  the  crystal  plane  and  the 
neutron  beam  will  be  designated  0.  The  interplanar  spacing  will  be 
represented  by  d  (see  Fig.  3).  Consider  two  parallel  waves  incident 
upon  the  crystal  planes  at  angle  8.  The  condition  necessary  for  con- 


Crystal 
Planes 


Fig.  3.  Waves  diffracted  from  crystal  planes 


structive  interference  to  occur  after  diffraction  is  given  by 


nA  =  AB  +  BC 


(3) 


where  n  is  an  integer  (1 ,2,3,...,n).  From  geometrical  considerations 
(Fig.  3) ,  it  is  seen  that 


AB  =  BC  =  d  sine. 


(4) 


Substituting  equation  4  into  equation  3,  the  familiar  Bragg  diffraction 
law,  equation  5,  is  obtained: 


ru  =  2d  sine.  (5) 

Scattering  which  obeys  equation  5  is  termed  "Bragg  scattering". 
2.3  Results  of  Different  Types  of  Scattering 

At  this  point,  mention  should  be  made  of  the  various  types  of 
scattering  which  will  be  of  interest.  All  of  the  available  information 
about  the  spectrum  will  be  obtained  from  the  coherent  elastic  scattering. 
Unfortunately,  incoherent  elastic  scattering  and  inelastic  scattering 
will  also  occur.  For  incoherent  elastic  scattering,  it  is  v/ell  known 
that  small  angle  scattering  is  much  more  probable  than  large  angle  scat- 
tering. On  the  other  hand,  inelastic  scattering  is  nearly  isotropic  for 
neutrons  having  an  energy  around  one  eV  or  below.  In  this  energy  range, 
most  of  the  commonly  used  diffraction  crystals  have  a  relatively  low 
absorption  cross  section.  Thus,  the  intensity  of  this  component  of  the 
background  will  be  roughly  proportional  to  the  volume  of  the  crystal 
which  is  intercepted  by  the  neutron  beam.  Naturally,  this  volume  will 
increase  as  the  difvraction  angle  decreases  up  to  the  point  where  the 
entire  crystal  is  in  the  neutron  beam.  From  this  point  on,  absorption 
may  cause  a  slight  decrease. 

As  a  result  of  these  two  effects,  one  may  expect  the  background  io 
be  a  function  both  of  angle  and  of  the  incident  neutron  spectrum.  For 
large  diffraction  angles,  the  background  should  be  relatively  constant. 
For  small  diffraction  angles,  however,  the  background  may  be  expected  to 
exhibit  a  strong  angular  dependence. 


One  further  comment  should  be  made  about  coherent  elastic  scattering, 
In  this  type  of  interaction,  the  neutron  is  interacting  with  an  atom 
which  is  tightly  bound  in  a  crystal  lattice.  Therefore,  the  neutron 
effectively  interacts  with  the  entire  crystal.  As  a  result  of  the  mass 
ratio  between  the  crystal  and  the  neutron,  the  neutron  loses  virtually 
no  energy  in  the  process  of  diffraction. 

2.4  Collimator  Effects  on  Resolution 

The  assumption  of  a  monoenergetic  parallel  beam  of  neutrons,  used 
in  deriving  equation  5 9  is  not  really  justified  in  neutron  diffraction 
work.  Instead,  the  neutron  beam  is  usually  colli mated  to  a  divergence 
of  somewhat  less  than  one  degree.  Also,  the  velocity  distribution  of  a 
thermal  neutron  emergent  from  a  nuclear  reactor  is  typically  (9,  10) 
described  by  a  modified  Maxwell i an  distribution  function  joined  to  a 
1/E  tail. 


in 
Intensity 


1/E  tail 


in   v 


Fig.  4.     Typical   thermal   leakage  spectrum 
from  a  nuclear  reactor 
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Since  measurements  will  be  taken  of  a  leakage  spectrum,  the  neutron 
current  density,  j(v),  in  the  emergent  beam  moy  be  represented  by 


j(v)  =/ 


3  ^V 

K.v^c       for  v<v 
1  c 


K0v~  for  v>v„ 
2         c 


(6) 


where  K.  and  K2  are  constants,  v  is  the  most  probable  neutron  velocity 
(in  the  sense  of  a  Maxwellian  distribution),  and  v  is  the  point  at  which 
the  two  functions  are  joined  to  best  fit  the  experimental  data. 

to  make  a  correction  for  the  finite  divergence  of  the  incident  neutron 
beam,  consider  Fig.  5.  Neutrons  from  a  point  source,  P,  may  strike  the 
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A  Y////////////?////////, 


Diffracting 
Crystal 


Fig.  5.  Effect  of  finite  divergence 
of  incident  neutron  beam 

crystal  surface  at  any  angle  between  a*  and  6".  If  L»W,  then 


60  =  (a*  -  a")  =  W/L. 


(7) 


However,  the  actual  neutron  source  at  the  entrance  to  the  collimator  is 
better  described  as  a  plane  source.  The  total  spread,  6eT,  of  the  inci- 
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dent  angle  of  the  neutron  beam  upon  the  crystal  is  then  given  by 

68T  =  f-   .  (8) 

This  may  easily  be  seen  by  roll  owing  point  P  as  it  travels  from  point 
A  to  point  B  across  the  end   of  the  collimator. 

To  convert  equation  8  to  a  more  suitable  form,  differentiate 
equation  5  to  obtain 

r>6\   =  2d  cose  60  .  (9) 

Equation  10  is  then  obtained  by  substituting  equation  8  into  equation  9: 

&x   =^cose  .  (10) 

Since  most  of  the  work  done  in  this  paper  is  in  terms  of  velocity  or  energy, 
equation  1  is  differentiated  and  substituted  into  equation  10  with  the 
following  result: 

2 
v  =   nhL  CCS  *  (11) 

Equation  11  yields  the  velocity  spread  introduced  into  the  diffracted 

beam  by  the  finite  divergence  of  the  incident  neutron  beam.  It  should  be 

noticed  that  the  velocity  spread  in  equation  11  is  net  constant  but,  indeed, 
depends  upon  the  angle  of  diffraction. 

2.5  Crystal  Effects 

2.5.1  Mosaic  Spread 

The  next  element  of  the  system  to  be  considered  is  the  diffraction 


1? 


crystal,  In  the  ,'.    don  of  equation  5,  a     ct  crystal  was  assumei  . 
Unfortun  tely,  all  large  single  crystals  are  com     of  a  large  number  of 
small  perfect  crystals.  These  are  usually  termed  crystallites.  Naturally, 
to  produce  a  large  single  crystal,  these  crystallites  must  have  their  cor- 
responding crystal  planes  highly  oriented.  The  degree  of  their  orientation 
is  usually  considered  (11)  in  the  following  manner.  Erect  a  line  perpendic- 
ular to  a  given  crystal  plane  in  each  crystallite  and  find  the  deviation  of 
each  of  these  normals  from  the  average  normal.  Plot  this  deviation  against 
its  frequency  of  occurrence.  A  roughly  Gaussian  curve  is  thereby  obtained. 
The  FWHM  of  this  curve  is  given  the  name  "mosaic  spread".  As  a  result  of 
this  type  of  crystalline  structure,  diffraction  will  occur  over  a  narrow 
band  of  angles,  even  with  a  perfectly  parallel  monochromatic  incident  neutron 
beam.  This  allows  a  curve  of  the  diffracted  beam  intensity  versus  diffrac- 
tion angle  to  be  cbtained.  Tne  area  underneath  this  curve  is  designated 
the  "integrated  reflectivity".  It  is  a  measure  of  the  number  of  neutrons 
having  a  given  energy  which  are  diffracted  from  the  crystal.  Bacon  (12) 
gives  the  following  relation  for  the  integrated  reflectivity,  R  : 

00 

r 


R  = 


aaA 


(1+a)  +  (l+2a)l/2  coth  A(l+2a)1/2 


(12) 


where  A  =  yt  /sine , 
r,e-A2/2n2 

yn»/2-rr 

Q  e  X3i!c2F2/sin26, 

y  is  the  macroscopic  aosorpticn  cross  section  for  tne  crystal, 

t  is  the  crystal  thickness, 
o 
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F  is  the  structure  factor  of   the  unit  cell, 

N  is  the  number  of  unit  cells  per  unit  volume,  and 

n  is  the  mosaic  spread  of  the  crystal. 
This  equation  is  valid  for  the  reflection  case  of  an  imperfect  crystal 
having  finite  absorption.  In  order  to  evaluate  equation  12  numerically,  the 
following  transformations  have  been  made:  let 


then  equation  12  becomes 


R  =  2/2n 


0 


Let 


then  by  rearrangement, 


1 i  = 


n/2 


6  =  A(l+2a)l/2 


V  = 


un/2TT 


Ve~u     du 


2  2 

(1+Ve"u  )  +  (l+2Ve"u  )'/2  cothG, 


T  =  H2Ve"u   , 


(13) 


2V. 


U    =    1/2    In    (~r) 


Equation  13  then  becomes 

1 
_  r 
V2n 


R  = 


1+2  V 


{1/2  Jln(|lr)}{(I+l)}+  T1/2  cothGs 


(14) 


14 


■1/2 
2/2"  n 


But,  G  =  AT  "  ,  leaving 

1+2V 
/ 

,6 


172 


(e^AT   -1)  dT 


T7T 


-.  (15) 


l 


(T+l)  (e2AT   -1)  tnCf^)  +4T1/2(e2AT   +1) 


With  benefit  of  hindsight,  the  following  relation  will  be  stated: 


TneBf<e>- 


(16) 


where 


f(e)  =  D  +  He1 


(16a) 


and  D,  H,  and  k  are  empirically  determined  constants. 

In  the  preceding  discussions,  it  has  been  shown  that  the  diffracted 
neutron  beam  will  not  be  monoenergetic,  but  will  instead  consist  of  a  range 
of  energies,  5E,  about  some  E.  The  practical  result  of  all  of  this  will  be 
to  raise  the  diffracted  beam  intensity,  but  to  reduce  the  beam  purity. 

It  has  been  mentioned  previously  that  the  intensity  and  the  divergence 
of  the  diffracted  beam  will  depend  upon  two  things.  These  are  the  mosaic 
spread  of  the  diffracting  crystal  and  the  divergence  of  the  incident  neutron 
beam.  Egelstaff  (13)  gives  the  following  equation  for  the  angular  width 
of  a  diffracted  beam: 


AG  = 


n2  ^  SQj.2 


1/2 


(17) 


He  also  suggests  that  for  a  diffracted  beam  of  given  divergence,  the  maxi- 
mum intensity  may  be  obtained  by  adjusting  the  collimation  of  the  neutron 
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beam  until  its  total  angular  spread,  66-r.  matches  the  mosaic  spread  of 
the  crystal. 

2.5.2  Beam  Interception 

The  last  important  characteristic  of  the  crystal  interaction  is 
entirely  one  of  geometry.  As  a  result  of  the  dimensions  of  the  incident 
beam  and  the  diffraction  crystal,  it  is  likely  that  for  small  angles  of 
incidence,  the  crystal  will  not  intercept  the  entire  beam.  Therefore,  a 
geometric  interception  factor,  G,  will- be  defined  as  follows: 


G  = 


1  if  sine  >  W/S 

,S  sine,  '.,     .       <  H,<  (18) 

( — g )  if  sine  <  W/S, 


where  S  is  the  width  of  the  diffraction  crystal. 

2.6  Neutron  Detector 

The  final  element  of  the  system  to  be  considered  is  the  neutron 
detector.  The  most  common  detector  used  in  this  work  is  a  large,  high 
pressure  BF~  probe.  Its  efficiency,  Eff,  may  be  rather  easily  calculated  (14) 
by  assuming  that  if  a  neutron  is  absorbed  in  the  detecting  medium,  a  signal 
will  be  obtained.  The  resulting  equation  is  simply 

Eff  =  1  -  exp(-zat),  (19) 

where  z.  is  the  macroscopic  absorption  cross  section  of  the  detecting 
a 

medium  and  t  is  the  thickness  of  the  detecting  medium  through  which  the 
neutrons  may  pass. 
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2.7  Theoretical  Correction  Factor 

Now  that  all  of  the  correction  factors  have  been  determined,  we 
are  able  to  calculate  a  theoretical  correction  factor,  C(v),  for  a  neutron 
spectrometer  where  C(v)  is  defined  by  equation  20: 

j(v)  =  C(v)  nQ(v),  (20) 

where  n  (v)   is  the  observed  countrate  at  any  given  setting  of  the  instru- 
ment, and  j(v)  is  the  corresponding  experimental  neutron  current.     C(v)  is 
given  by  the  proper  combination  of  equations  11,  16,  18,  and  19.     Thus,  we 
have 

c<v>  "  2  Wvncose    f<6>  «'  "  e"Eat>  <21> 

2.8    Cross  Section  Determination 

The  determination  of  a  neutron  cross  section  is,  in  theory,  a 

relatively  simple  matter.     When  a  sample  is  placed  in  a  neutron  beam,  (15,16) 

the  new  beam  intensity,  nT,  is  given  by 

-Not 
nT  =  nQe    n     s,  (22) 

where  n„  is  the  unattenuated  countrate,  N  is  the  number  of  nuclei  per 
o  n 

unit  volume  of  the  sample,  a  is  the  microscopic  neutron  cross  section,  and 
t  is  the  thickness  of  the  sample.  Equation  22  may  be  written  in  the  some- 


what more  convenient  form 


0  ■  o-*"^  •  (»> 

n  s    T 
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Thus,  the  velocity  dependence  of  a  may  be  readily  determined.  Naturally, 
the  particular  cross  section  laeasured  depends  upon  the  geometry  of  the 
experiment.  Typically,  the  total  cross  section  will  be  measured  in  a 
"good  geometry"  experiment. 

With  the  exception  of  specially  prepared  specimens,  most  crystalline 
materials  consist  of  a  great  number  of  randomly  oriented  microcrystals. 
If  a  suitable  sample  is  placed  in  a  collimated  neutron  beam  having  a 
continuous  spectrum,  diffraction  may  occur  from  each  microcrystal  accord- 
ing to  equation  5.  Equation  5,  however,  predicts  a  maximum  wavelength, 
which  may  be  scattered  for  any  given  set  of  crystal  planes.  Under  the 
limiting  conditions  of  maximum  wavelength,  equation  5  becomes 


*  =  2d 
c 


(24) 


and  is  sometimes  known  as  the  Bragg  cutoff  formula. 

Thus,  it  is  apparent  that  neutrons  having  a  shorter  wavelength  than 
given  in  equation  24  will  undergo  more  scattering  than  those  having 
longer  wavelengths.  Hence,  one  expects  a  sharp  change  in  the  scattering 
cross  section  at  certain  energies.  This  effect  is  very   pronounced  in  a 
number  of  materials,  as  is  shown  in  Fig.  6. 
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Fig.  6.  Bragg  cutoff  wavelengths  of 
different  materials  (17) 


3.0  EXPERIMENTAL 
3.1  Description  of  the  Reactor 

The  Kansas  State  University  TRI6A  Mark  II  (KSUTMII)  nuclear  reactor 
is  a  light  water  cooled,  graphite  reflected  reactor.  The  fuel  consists 
of  a  uranium- zirconium  hydride  alloy  with  a  20%  enrichment  in  U-235. 
Because  of  the  nature  of  the   hydrogen  bonding  in  the  fuel,  nearly  all  of 
the  moderation  for  neutron  energies  below  .13  eV  occurs  in  the  cooling 
water  surrounding  the  fuel  elements.  Above  that  energy,  moderation  may 

.3  place  in  either  the  fuel  or  the  cooling  water.  At  the  present  time, 
the  maximum  licensed  power  level  is  100  kw.  The  varied  experimental 
facilities  are  indicated  in  Fig.  7.  They  include  four  beamports  which 
are  designated  by   the  direction  in  which  they  face,  i.e.,  S.E.,  N.E.,  N.W., 
and  S.W. 

The  S.E.  beamport  looks  at  the  graphite  reflector  tangentially  to 
2  reactor  core.  Therefore,  all  neutrons  and  core  gammas  emergent  from 
this  port  must  have  been  scattered  at  least  once.  After  this  first  scat- 
tering, most  of  them  must  still  pass  through  a  large  thickness  of  graphite 
before  entering  the  beam  tube.  Thus,  this  port  is  used  for  experiments  requir- 
ing a  well -moderated  beam  of  neutrons  and  a  high  neutron  to  gamma  ratio. 

The  N.E.  beamport  looks  directly  at  the  outside  of  the  reactor  core 
through  a  hole  in  the  graphite  reflector.  This  port,  then,  will  yield  a 
relatively  high  flux  of  epi thermal  and  fast  neutrons  together  with  a  high 
gamma  dose  rate. 

The  N.W.  and  S.W.  ports  are  identical  except  for  their  location.  Like 
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the  N.E.  port,  they  look  radially  at  the  outside  of  the  reflector.  However, 
they  are  not  aligned  with  any  hole  in  the  reflector.  Therefore,  all  neutrons 
emerging  from  these  ports  must  have  passed  through  at  least  12  inches  of 
graphite. 

Table  I  gives  the  best  available  estimate  of  the  relative  neutron  fluxes 
and  gamma  dose  rates  at  selected  points  in  the  reactor. 

Table  I.  Radiation  Fields  at  Selected  Points 

in  Reactor  for  100  kw  Power  Level  (18) 


Position 


Neutron 
(n/cm  -sec) 


Gamma 
(rads/sec) 


Fast  (>10  KeV)   Thermal  (<.21  eV) 


N.E.  port  next       ■  8.0  x  10 
to  core 

Reactor  pool  outside    2.7  x  10 
reflector 


11 


10 


8.0  x  10 


2.7  x  10 


11 


11 


4.0  x  10' 


1.8  x  10' 


All  of  the  beam  tubes  terminate  at  the  outer  edge  of  the  reflector. 
Because  of  manufacturing  and  constructional  tolerances,  there  is  approximately 
1/2  inch  clearance  between  the  ends  of  the  beam  tubes  and  the  outer  surface  of 
the  reflector.  The  space  in  between  these  two  surfaces  is  filled  with  alum- 
inum shims  which  are  curved  on  one  side  to  match  the  curvature  of  the  reflector 
housing.  Because  the  reflector  is  canned  in  1/2  inch  thick  aluminum  and  the 
ends  of  the  beam  tubes  are  about  5/32  inch  thick,  neutrons  must  penetrate  a 
substantial  amount  of  aluminum  before  they  enter  the  beam  tube.  Since  the  mean 
free  path  for  thermal  neutrons  in  aluminum  is  about  two  inches,  one  may  expect 
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a  certain  perturbation  in  the  thermal  neutron  spectrum  corresponding  to 
the  variation  in  the  aluminum  cross  section.  The  exact  magnitude  of  the 
perturbation  is,  however,  wery   difficult  to  predict. 

Fig.  8  shows  the  standard  fuel  loading.  This  particular  "loading  is 
designated  number  6.  Note  that  there  are  four  fuel  elements  in  the  F  ring. 
These  fuel  elements  may  be  placed  in  any  convenient  location  in  this  ring. 
Also,  up  to  six  of  the  graphite  reflector  elements  may  be  removed  and  still. 
maintain  enough  reactivity  to  attain  a  power  level  of  100  kw,  provided  that 
no  extraneous  poisons  are  inserted  into  the  core.  Figures  9,  10,  and  11 
show  fuel  loadings  which  have  been  modified  to  achieve  certain  desired 
results.  All  comparisons  made  in  the  following  descriptions  of  the  modified 
fuel  loadings  will  be  with  respect  to  the  standard  fuel  loading. 

In  Fig.  9,  the  four  fuel  elements  in  the  F  ring  have  been  placed  in 
a  position  close  to  the  tangential  port  in  order  to  increase  the  neutron 
flux  from  this  port.  This  loading  is  designated  6'. 

Fig.  10  indicates  a  fuel  loading ,  6"  which  will  increase  the  fast  and 
epi thermal  nautron  fluxes  from  the  N.E.  port.  Note  that  this  loading  permits 
a  minimum  of  water  between  the  fuel  and  the  reflector  housing.  This  is  yary 
important  because  of  the  moderating  properties  of  water. 

Fig.  11  indicates  a  fuel  loading,  6,f,9  which  will  increase  the  thermal- 
ization  of  the  neutron  beam  from  the  N.E.  port.  Note  that  the  essential 
change  is  that  'the  graphite  reflector  elements  have  been  replaced  by  a  much 
more  effective  moderator,  water. 


*The  designation  given  to  the  different  fuel  loadings  is  that  given  by  the 
reactor  supervisor  and  entered  into  the  reactor  log  book. 
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Fig.  8.  6  core  load 
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N.E.  Port 


■  Core  Loading  Diagram 
■  F   Fuel  Element 
TC   Instrumented  Fuel  Element 
G   Graphite  Reflector  Element 

Fig.  9.  6'  core  load 
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Fig.  10.  6"  core  load 
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Fig.  11.  6"'  core  load 
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3.2  Description  of  the  Collimator 

The  radiation  beam  emergent  from  a  completely  open  beam  tube  is 
much  too  hazardous  and  unwieldy  for  most  experiments.  As  a  result, 
collimators  are  used  to  reduce  the  radiation  beam  to  more  manageable 
dimensions  and  to  provide  a  means  of  reducing  its  angular  divergence. 
The  collimator  must  therefore  fill  the  beam  tube  for  a  distance  of 
several  feet.  The  only  opening*  will  be  a  suitably  sized  hole  through 
which  a  radiation  beam  of  the  desired  cross  section  is  allowed  to  pass. 
Because  of  the  sizable  gamma  and  fast  neutron  dose  rates  present  near 
reactor  cores,  most  collimators  contain  sizable  quantities  of  lead,  steel 
and  concrete.  As  a  result,  the  typical  collimator  is  a  rather  large, 
heavy  device. 

Since  no  collimator  was  available  for  the  KSUTMII,  one  had  to  be 
constructed.  The  design  requirements  were  as  follows: 

1.  The  collimator  must  fit  in  all  four  beam  tubes  and  allow 
the  safety  shutter  to  be  completely  closed  and  the  outer 
door  to  be  sealed  and  locked, 

2.  Shielding  plugs  must  be  provided  so  that  the  reactor  may 
be  run  at  100  kw  with  the  collimator  in  any  beam  tube  and 
still  not  produce  a  radiation  hazard  when  the  collimator 
is  not  in  use, 


*0ther  openings  for  instrumentation,  etc.,  are  sometimes  present,  but  they 
are  usually  stepped  in  some  manner  to  avoid  a  radiation  hazard. 
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3.  The  collimator  must  provide  good  collimation  for  fast  neutrons, 
thermal  neutrons,  and  gamma  rays, 

4.  The  collimator  must  not  present  a  storage  problem  because 
of  neutron  activation  of  materials  having  long  half  lives, 

5.  The  collimator  must  provide  an  air  seal  when  in  place  so  as  to 
prevent  hazardous  quantities  of  radioactive  gases  from  venting 
from  the  beam  tube  into  the  reactor  bay,  and 

6.  The  collimator  must  be  constructed  from  inexpensive  and  avail- 
able materials  so  far  as  possible. 

In  order  to  determine  a  size  for  the  proposed  collimator,  it  was  first 
necessary  to  check  the  dimensions  of  all  beamports.  The  diameters  were 
found  in  an  appropriate  publication  (19),  but  no  reference  could  be  found 
which  gave  the  linear,  dimensions.  Thus,  all  of  the  beam  tube  plugs  had  to 
be  removed  and  the  beam  tubes  measured.* 

The  internal  dimensions  of  the  hole  through  the  collimator  were  set  by 
the  size  of  the  diffraction  crystal  and  the  size  of  the  available  lining 
tube. 

At  that  time,  the  maximum  power  level  for  which  experimental  dose  rate 
measurements  from  an  open  beamport  were  available  was  three  watts.  Thus, 
a  series  of  measurements*  was  prepared  for  different  power  levels  at  the 
S.E.,  N.E.,  and  N.W.  ports.  The  N.W.  port  was  chosen  in  preference  to  the 
S.W.  port  because  of  the  arrangement  of  equipment  within  the  reactor  bay. 
After  these  data  had  been  evaluated  and  more  confidence  gained  in  the  meas- 
uring techniques,  the  N.E.  port  was  completely  opened  and  the  reactor  was 


*See  Appendix  A  for  a  summary  of  these  measurements. 


operated  at  higher  power  levels  to  determine  the  magnitude*  of  the  scattered 
radiation  at  various  points  around  the  reactor  bay. 

The  aforementioned  series  of  measurements  yielded  enough  data  to  enable 
the  collimator  to  be  designed  and  constructed.  Fig.  12  shows  the  internal 
construction  of  the  collimator  while  Fig.  13  is  a  photograph  of  the  actual 
coll"'  lator  and  its  associated  equipment. 

To  simplify  geometrical  considerations,  it  was  decided  to  construct 
the  collimator  to  produce  a  radiation  beam  having  a  rectangular  cross  sec- 
tion. A  suitable  tube  and  sizable  quantities  of  lead  and  polyethylene  were 
obtained  from  government  surplus.  Simple  calculations  indicated  that  nearly 
10  inches  of  lead  would  be  required  to  reduce  core  gammas  and  capture  gammas 
to  the  desired  dose  rate.  The  lead  was  cast  into  five  two-inch  discs  in 
order  to  better  distribute  the  weight  and  to  conform  to  the  "removal  cross 
section"  concept.  The  lead  discs  were  separated  by  borated  paraffin  which 
was  cast  in  place. 

Cadmium  and  polyethylene  discs  were  placed  on  the  inside  end  of  the 
collimator  to  thermalize  intermediate  and  some  fast  neutrons  and  to  prevent 
the  thermal  neutrons  from  activating  the  impurities  in  the  lead.  Although 
these  impurities  may  become  highly  radioactive""  under  neutron  bombardment, 
they  appear  to  be  nearly  pure  beta  emitters  with  a  fairly  short  half  life. 
Because  of  the  low  neutron  flux  in  which  the  collimator  is  placed  and  the 
amount  of  self-shielding  provided  in  the  collimator  design,  it  was  decided 


'See  Appendix  A  for  a  summary  of  these  measurements. 
'See  Appendix  B. 
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that  these  impurities  would  cause  no  storage  problems.  Also,  an  additional 
amount  of  boric  acid  was  placed  against  each  face  of  the  lead  discs  during 
construction  of  the  collimator  as  an  added  precaution. 

As  is  obvious  from  Figures  12  and  13,  the  collimator  consists  of  two 
major  parts.  The  larger  part  fits  into  the  beam  tube  and  may  contain  the 
shielding  plugs.  The  smaller  part  is  an  extension  to  provide  better  col- 
limation  for  thermal  neutrons.  It  also  serves  as  a  handle  with  which  to 
insert  and  to  extract  the  collimator  from  the  beam  tubes.  When  this  exten- 
sion (snout)  is'  removed,  the  beam  port  safety  shutter  and  outer  door  may  be 
completely  closed  with  the  main  body  of  the  collimator  remaining  in  the  beam 
tube.  The  shielding  plugs  fit  tightly  enough  to  serve  the  same  purpose  if 
the  snout  is  removed. 

One  special  precaution  should  be  observed  when  moving  the  collimator. 
The  square  tube  along  the  longitudinal  axis  of  the  collimator  is  a  magnesium- 
aluminum  alloy.  All  welded  joints  along  this  piece  are  fragile  and  should 
be  bumped  as  little  as  possible.  Special  care  should  be  taken  when  insert- 
ing the  collimator  into  a  beam  tube. 

Three  tubes  were  run  through  the  length  of  the  collimator  to  allow 
future  expansion.  They  may  be  used  for  such  things  as  instrumentation  cables 
and  refrigerants. 

3.3  Description  of  the  Neutron  Spectrometer 

To  insure  adequate  counting  statistics,  the  cross-sectional  area  of 
the  neutron  beam  being  diffracted  is  fairly  large,  typically  on  the  order 


of  a  few  square  inches.  In  order  for  the  detector  to  intercept  the  dif- 
fracted beam,  it  usually  has  a  cross  section  of  at  least  the  same  dimen- 
sions. Since  BFn  detectors  are  generally  used  because  of  their  ease  of 
0]  .ration  and  their  excellent  gamma  rejection,  the  detector  must  be 
nearly  20  inches  long  in  order  to  achieve  the  desired  efficiency  for  end- 
on  counting.  Because  of  its  large  size,  the  detector  must  be  efficiently 
~lded  to  eliminate  any  room  scattering  contributions  to  the  background. 
This  shield,  then,  will  ..eigh  several  hundred  pounds. 

To  obtain  good  angular  resolution  with  such  a  large  neutron  beam, 
the  detector  is  often  placed  some  distance  from  the  diffractin  crystal. 

:,  a  good  neutron  spectrometer  must  be  capable  of  accurately  supporting 
a  long  arm  with  a  heavy  weight  on  the  far  end.  This  is  the  reason  for  the 
large  size  of  such  instruments. 

The  design  of  the  KSU  neutron  spectrometer  was  determined  by  the 
following  considerations: 

1.  Available  neutron  flux  at  the  exit  of  the  collimator, 

d. .   oOS  C  -j 

3.  Usable  at  the  S.E.,  N.E.,  and  N.W.  beam  ports, 

4.  Accuracy, 

5.  Ease  of  operation, 

6.  Energy  range  covered  and  purity  of  the  diffracted  beam,  and 

7.  Further  possibility  of  expansion  of  the  capabilities  of  the  instrument. 
Preliminary  measurements  of  the  neutron  flux  emergent  from  the  beamport 

collimator  at  the  S.E.  port  yielded  a  value  for  the  neutron  flux  at  the 

2 

position  of  the  detector  of  10  n/cm  -sec/watt.  It  was  thus  concluded  that 
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a  neutron  spectrometer  could  be  constructed  which  could  give  meaningful 
measurements. 

Exclusive  of  labor  costs,  the   neutron  spectrometer  and  the  collimator 
were  constructed  for  approximately  $1500.  Much  use  was  made  of  surplus 
materials  but  such  things  as  precision  gears  and  diffracting  crystals  had 
to  be  obtained  through  normal  commercial  channels. 

The  maximum  length  of  the  detector  arm  was  determined  by  the  space 
between  the  reactor  shielding  and  the  shielding  wall  for  the  neutron  gen- 
erator located  near  the  N„E.  beam  port.  This  is  a  distance  of  12  feet. 
The  height  of  the  diffracting  crystal  was  determined  by  the  distance  of 
the  center  of  the  beam  ports  above  the  reactor  bay  floor.  This  was  a  major 
design  consideration  as  the  distance  is  only  36  inches. 

Accuracy  proved  to  be  the  major  problem.  Although  the  detector  could 
be  positioned  within  1/20  of  a  degree  of  the   desired  setting3  reproducib'.a 
results  could  not  be  obtained.  Early  measurements  indicated  that  movement 
of  the  arm  produced  enough  internal  strain  to  tilt  the  crystal.  The  main 
support  ring  was  investigated  and  found  to  have  a  saddle-shaped  surface. 
This  surface  was  then  finished  to  a  tolerance  of  ±.002".  The  main  bearings 
and  the  gearbox  mountings  were  redesigned.  These  s\;eps  helped  to  correct 
the  problems  but  did  not  entirely  eliminate  them.  At  this  time  a  mere  sen- 
sitive optical  method  of  checking  the  internal  alignment  was  devised. 

An  optical  level  was  placed  about  15  feet  from  the  crystal  holder.  A 
small  light  bulb  was  fastened  on  the  outside  end  of  the  detector  shield  and 
a  front-surfaced  mirror  placed  in  the  crystal  mounting.  The  crystal  table 
was  rotated  until  the  light  was  observed  to  be  centered  on  the  cross  hair 
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scratched  on  the  front  surface  of  the  mirror.  The  arm  was  then  rotated 
through  an  angle  of  approximately  100  degrees  (the  operating  range  of  the 
instrument)  and  the  apparent  displacement  of  the  light  note  . 

This  optical  method  allowed  the  causes  of  the  irregular  behavior  of 
the  crystal  to  be  finally  determined.  They  were  found  to  be  a  movement  of 
the  crystal  table  caused  by  some  electrical  leads,  and  deflections  of  the 
support  ring  under  the  weight  of  the  main  arm  and  its  associated  counter- 
weights. The  use  of  more  flexible  wires  and  the  placement  of  hydraulic  jacks 
at  strategic  points  under  the  support  ring  cured  these  remaining  difficulties. 

The  gearbox  provides  the  necessary  2:1  reduction  in  angular  displacement 
between  the  detector  and  the  diffracting  crystal.  A  photograph  of  the 
gearbox  and  the  crystal  table  assembly  is  shown  in  Fig.  16.  The  gearbox 
is  fastened  to  the  upper  main  arm  bearing  by  three  10-24  cap  screws.  These 
screws  are  easily  visible  on  the  bottom  shelf  of  the  gearbox.  Before  removing 
these  screws,  it  is  first  necessary  to  loosen  the  bottom  left  gear.  After 
doing  this,  the  three  screws  may  be  removed.  The  entire  gearbox  and  crystal 
table  assembly  may  then  be  removed  by  lifting  them  straight  up.  If  only  the 
crystal  table  assembly  is  to  be  removed,  it  is  necessary  to  loosen  only  the 
gear  in  the  upper  left  of  the  gearbox.  The  entire  crystal  table  assembly 
may  then  be  removed  by  lifting  it  straight  up. 

The  manner  in  which  the  gearbox  achieves  its  2:1  reduction  is  not  entirely 
obvious  from  Fig.  16.  Although  the  gearbox  itself  is  fastened  to  the  upper 
main  arm  bearing,  the  small  gear  in  the  lower  left  of  the  gearbox  is  attached 
to  a  pin.  This  pin  is  press  fitted  into  the  main  shaft  which,  in  turn,  is 
firmly  attached  to  the  base  of  the  instrument.  Thus,  this  gear  remains  sta- 
tionary and  the  entire  gearbox  rotates  around  it.  The  pin  ends  just  below 
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Fig.  16.  Gearbox  and  crystal  table  assembly 


-  middle  shelf  in  the  gearbox.  The  shaft  supporting  the  crystal  table 
assembly  then  rests  on  the  top  of  the  pin  for  vertical  support.  The  shaft 

the  right  side  of  the  gearbox  is  solid  and  only  transmits  the  motion 
between  the  lower  and  upper  sets  of  gears.  Since  the  ratio  of  the  t\ 
gears  in  the  top  is  1:1,  the  desired  overall  ratio  is  achieved. 

3ec^-^c  of  the  type  of  experiments  which  were  to  be        !,  as  well 
as  constructional  problems,  it  was  decided  not  to  construct  a  massive 
shield  arou.-.c  -che  diffraction  crystal.  Instead,  an  external  beam  catcher 
to  be  used.  Also,  it  was  felt  that  this  approach  would  help  to  reduce 
e  background  for  small  diffraction  angles.  Later  measurements  of  both 
.  scattered  radiation  within  the  reactor  bay  and  of  the  collimated  be 
licated  that  no  particular  radiation  safety  problems  would  be  presented 
by  this  decision. 

It  was  felt  that  automation  of  the  entire  instrument  would  be  too 
costly  and  add  greatly  to  the  undertaking.  Therefore,  only  the     -  lent 
which  provided  rotation  of  the  crystal  about  the  vertical  axis  was  eqi 
with  a  motor  drive.  This  was  done  mainly  because  of  the  fineness  of  the 
adjustment.  The  reversible  motor  is  manually  controlled  with  a  three- 
positions  center-off,  switch.  The  motor  is  coupled  to  the  adjustment  thumb- 
screw with  a  suitably  sized  0-ring. 

Two  other  adjustments  are  provided  for  the  crystal.  They  allow  ro- 
tation abouj  a  horizontal  axis  as  well  as  translation  of  the  entire  crystal, 
Complete  instructions  for  the  alignment  of  he  neutron  spectrometer  are 
provided  in  the  following  section. 

In  order  to  minimize  the  effects  of  fluctuation  in  the  reactor  power, 
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a  monitor  assembly  was  devised.  Initial  attempts  at  placing  a  partially 
shielded  BF0  probe  near  the  collimator  met  with  failure.  The  monitor 
probe  assembly  so  constructed  proved  to  be  highly  dependent  upon  the 
main  arm  position.  The  BF~  probe  was  finally  placed  within  a  polyethylene 
block  which  was  placed  between  two  masonite  spacers  (see  Fig.  17).  The 
outer  wooden  shielding  plug  was  removed  from  another  beam  tube  and  the 
monitor  assembly  was  then  inserted  to  a  depth  of  14  1/2  inches.  Two  addi- 
tional polyethylene  blocks  were  placed  between  the  monitor  probe  and  the 
inner  beam  plug  to  help  thermal ize  some  of  the  fast  neutrons.  All  poly- 
ethylene blocks  used  were  cylinders  9  inches  long  by  5  3/4  inches  in 
diameter. 

The  electronic  equipment  used  was  connected  as  indicated  in  Fig.  18. 
The  idea  is  to  collect  a  preset  number  of  counts  from  the  diffracted  beam 
(data  channel).  When  this  preset  number  is  reached.,  the  timer  and  the 
power  monitor  (monitor  channel)  counter  are  stopped.  Thus,  the  monitor 
channel  yields  a  number  which  is  proportional  to  the  time  integral  of 
the  reactor  power  during  the  counting  time. 

The  timer  allows  one  to  make  a  resolving  time  correction  to  the 
data  channel  information.  No  such  correction  for  the  monitor  channel 
is  necessary  so  long  as  the  reactor  is  run  at  a  relatively  constant  power. 
By  suitable  adjustment  of  the  position  of  the  monitor  probe  within  the 
beam  tube,  the  countrate  of  the  monitor  channel  can  be  made  to  greatly 
exceed  that  of  the  data  channel.  (Care  must  be  taken  not  to  saturate 
the  monitor  channel  scaler,  however.)  Thus,  the  statistics  of  the  final 
data  will  be  determined  mainly  by  those  of  the  data  channel. 
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is  not  used  in  obtaining  the  final  data,  a     it  was 
logue  voltage  proportional  to  the  an    .  ._. 
The  -am  is  sha     rig.  19.  It  is  a  simple 

circuit  '  .   voltage  regulator.  This  circ  to   used 

conjunction  with  a  TMC  GAMMA  SCOPE  operating  in  the  Mossbam  r  mode. 


□  n,  h  watt.. 
AAA/ 


12  volts 


1  N3020B 


10  volt, 

5%) 


loco  a 


100  :. 

IRC  Type  151 


Fig.  19.  Analogue  position  circuit 


_-;.r  crystals  were  selected  for  the  diffraction  grating  because  of 
their  cost,       lity,  and  suitability.  The  maximum  size  available 
2  inches  by  1  inch  with  a  thickness  of  6  millimeters.  Two  crystals 
were  obtained  with  a  face  cut  parallel  to  the  (HI)  plane  and  two  cut 

rallel  to  the  (200)  plane.  The  (111)  crystals  have  the  advantage  (20) 
that  they  are  known  to  produce  no  second  order  diffraction  pattern.  They 
do,  however,  produce  a  sizable  third  order  pattern.  The  (200)  crystals  ■ 
have  (21)  a  second  order  diffraction  pattern,  but  they  also  have  a  smaller 
interplanar  spacing  and  a  smaller  third  order  diffraction  pattern  than  the 
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(111)  crystals.  Also,  the  (200)  crystals  have  the  advantage  of  a  smaller 
mosaic  spread  and  therefore ,  can  be  used  to  make  higher  resolution  meas- 
urements than  are  possible  with  the  (111)  crystals. 


3.  t  Procedure 

ajor  Adjustments 

Any  adjus       -  .;/  indicated  in 

ins       s  should  not      sction  of 

iced  _j-      Even  t ,  ...    justi  ent  of  the  central 

shaft  should  be  undertaken  only  after  great  deliberat'  .  .  If  it 
is  2nt,  it    t  be  reset        icial  jig  ai 

a  dial  indicator. 
II.  Align  ent  with  tl  \        ctor  Shut  Do 

.'is  procedure,  there  should  be  no  excuse  for  anyone 

to  receive  s ficant     tion  dose.  It  has  :>e    .  Dnstrated 

(see  t     etc  r   .    rt  log)  that  the  cose  rate  .       from  t 
beamport  may  be  as  high  as  10  ?/'r.r.   (a:r':er  a  10(    'un)f  even 
thoi     _    ;tor  is  subcritical.  Also,  there  may  be  a  sizable 

)unt  of  radioactive  z&s   generated  inside  the  fc    Drt.  Since  be 
of  these  hazards  rapidly  cecrease  with  time,  they  may  be  reduced  "co 
minor  proportions  by  mer  .   liting  the  proper  1     of  time. 

1.  Consult  with  the     ;    .  ervisor  about  the  ast  operating 
history  c      .  ctor  and  the      -  techniques  for  opening 
the  beam  ports. 

2.  Obtain  suitable  radiation  measuring  equipment,  an  automatic  level 

or  transit,  a  carpenter's  level,  a  light  sourc  ■.  

inserted  into  the  collimator,  and  the  special  front-surfaced 

irror. 
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3.  Using  suitable  technique,  remove  the  beamport  plugs  and  "insert 
the  collimator  into  the  chosen  beamport.  Make  the  proper  entries 
in  -  te  Lamport  log. 

4.  Place  the  spectrometer  assembly  roughly  in  position  with  the  coge 

of  the  crystal  table  about  1"  from  the  front  of  the  collimator  snout. 
Place  the  carpenter's  level  across  the  bottom  ring  and  adjust  the 
three  jackscrews  until  the  instrument  is  level.  Co  not  raise  the 
instrument  any  higher  than  necessary.   (See  Fig.  20) 

5.  Remove  the  snout  and  the  internal  shielding  plugs  from  the  collimator. 
Insert  the  light  source  until  it  hits  the  inside  end  of  the  collimator. 

6.  Using  the  automatic  level  or  transit,  move  the  whole  spectrometer 
assembly  until  center  of  the  crystal  holder  is  in  line  with  the 
collimator  (vertical  axis  only).  The  overhead  crane  will  be  necessary 
for  this  movement, 

7.  Place  1  drop  of  oil  on  each  bearing  in  the  gearbox.  Generously  oil 
both  main  shaft  bearings  on  the  main  arm. 

8.  Relevel  the  spectrometer  by   adjusting  the  three  jackscrews  and  acjust 
height  of  the  instrument  until  the  bottom  of  a  neutron  beam  from  the 
collimator  would  just  strike  the  top  of  the  crystal  cajle. 

9.  Insert  the  front-surfaced  mirror  into  the  crystal  holder. 

10.  Using  the  transit  or  automatic  level,  rotate  the  spectrometer  arm  and 
make  sure  that  the  center  of   rotation  of  the  crystal  holder  is  correct, 
The  upper  crystal  table  may  be  slid  back  and  forth  to  correct  for  any 
misalignment. 

11.  ;;ove  the  light  source  in  the  collimator  sufficiently  close  to  the 
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outside  so  that  a  usable  reflection  from  the  mirror  may  be  obtained. 
Place  the  arm  at  approximately  a  30°  angle  and  adjust  iha   crystal 
holder  ur.c: ':   the  spot  of  light  illuminates  the  detector. 

12.  Remove  the  mirror  from  the  crystal  holder  and  the  light  from  the 
collimator.  Move  the  beam  definition  shields  on  front  to  produce 
a  neutron  beam  h"   wide  and  Vs  high.  Make  sure  that  the  beam  will 
strike  on'ly   the  crystal  and  only  in  the  center.  Position  the  bea~. 

;cher  behind  the  spectrometer. 

13.  The  collimator  is  so  designed  so  as  to  enable  it  to  be  placed  in 
a.vy  of  the  four  beamports  and  have  both  the  shutter  and  the  outer 
door  closed.  If  this  is  done,  however,  the  collimator  should  have 
its  own  internal  shielding  plugs  in  place.  With  these  in  place, 
the  radiatic;.  field  outside  the  beamport  will  be  less  than  2  mr/hr. 
even  with  the  reactor  operating  at  100  kw.  Hence,  if  the  spectrometer 
is  to  be  left  unattended  for  any  length  of  time,  the  collimator 
should  Le  plugged  to  avoid  accidents.  Also,  remember  to  bring 

the  beamport  status  sign  up  to  date. 
III.   Neutron  Alignment 

1.  After  checking  with  the  Reactor  Supervisor.,  remove  the  shielding 
plugs  from  the  collimator  and  attach  the  collimator  snout.  Hang  the 
RAM  probe  over  the  snout  and  lock  it  in  position  with  the  beamport 

1  cck. 

2.  fter  allowing  a  suitable  warmup  time  for  the  associated  electronic 
equipment,  take  the  reactor  up  to  a  power  level  of 
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beampcrt  power  level         maximum  expected  countrate 

a.  tangential  100  watts  36,000  cpm 

b.  fast  30  watts  35,000  cpm 

c.  regular  100  watts  42,C0C  cpm 

3.  Without  a  crystal  in  place  move  the  spectrometer  arm  until   the 
detector  is  barely  out  of  the  direct  neutron  beam.     Take  a  one 
minute  count  and  move  the  arm  \  cm  closer  to  the  beam.     Continue 
until  the  width  of  the  beam  has  been  traversed.     By  plotting  these 
data,  the  zero  point  on  the  calibration  ring  may  be  determined. 

4.  Either. scram  the  reactor  or  reduce  the  power  level  to  one  watt 
or  less. 

5.  Carefully  install  the  proper  crystal  in  the  crystal  mount  and 
then  move  the  spectrometer  arm  until  it  makes  an  angle  of  about 
30°  with  the  emergent  neutron  beam. 

6.  Increase  the  reactor  power  to 

beampcrt 

a.  tangential 

b.  fast 

c.  regular 

7.  Being  wery   careful  not  to  get  any  portion  of  the  body  in  direct 

■  contact  with  the  radiation  beam,  use  one  hand  to  adjust  the  crystal 
using  the  adjusting  screws  underneath  the  crystal  table  until  the 
maximum  countrate  is  obtained.   (See  Fig.  21) 
3.  Rotate  the  crystal  about  a  vertical  axis  until  only  background  is 

obtained.  Rotate  1°  more  in  same  direction. 
9.  Increase  reactor  power  to  triple  that  indicated  in  step  6  above. 
10.  Using  a  ratemeter  with  a  time  constant  of  1  second  and  a  recorder, 
rotate  the  crystal  table  back  through  the  peak  countrate  position 


power 

level 

10 

kw 

3 

kw 

10 

kw 
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until  only  background  is  beinc        :orded.     Use  the  i   ■  -.ive 

iture  for  this  part.     The  speed  c?  the  ..  irive  may  be 

.armine^.  b;     ti  ling  the  :ings  c  .       .  edge  of 

ystal  table.     This  procedure  gives  the         .  ng  cur. 
a  measur_  of  the  instrument  resolution. 

11.  Reverse  t    -  direction  of  rotation  of  the  crystal   t     le  until 

;rate  again  reaches  the  peak  and  then  turn  off  the 
;or. 

12.  Position  an?,  for  first  reading. 

13. 

tangential  port        IOC 
fast  port 
c.  regular  port  100  kw 

and  .\    t  reactor  operator  to  swi     d  automatic  control. 

When  power  lava!  has  stabilized  (5  minutes)  de    ine  neutron 

specti     taking  a  certain  number  of  counts  for  each  of 

spectrometer  ar 

15.  Upon  completion  of  data  taking,  shut  down  t    :.c;;/~  and  allow 

it  to  cool  until  such  iime  as  it  is  safe  to  reinsert  the  shielding 

plugs  into  the  collimator. 

. . .  Radi ati on  Moni tori  ng 

.  .  ing  part  II  of  these  alignment  instructions,  the  only  monitoring 

necessary  should  ba  that  of  the  gamma  field  during  placement  of  the 

collimator.  If  proper  time  has  been  allowed  for  the  reactor  core  to 

cool9  then  the  radiation  field  should  be  at  most  40  mr/hr  gamma.  Fr 

the  time  and  th      :as  involved  it  would  seam  highly  unlikely  that 

an  exzzz^ra   of       n  40  mr  would  result.  Indeed,  1  mr  would  be 
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more  likely. 

.    rt  III  will  require  somewhat  more  monitoring.     II 1—1  tl 
III-6  should  onl\    require  an  observer  to  make  sure  that  the  person 
making  the  ad;  s     ys  out  of  the   ii  \        'a<    ation  beam. 

roperly  performed,  no  one  would  come  within  one  foot  of  the  direct 
radiation  beam.     Ill— 7  and  III— 8  will   require  that  a  finger  dosimeter 
be  worn.     The  radiation  field  underneath  the  crystal  table  is  about 
w  mr/hr  (at  100  kw)  and  results  mainly  from  the  Co  strips  on  the 
front  of  the  collimator.     Even  so9  it  would  be  unlikely  that  the  hands 
would  receive  more  than  20  rnr.  in  the  necessary  time.     If  the  hands 
are  placed  in  the  direct  radiation  beam,  however,  that  is       Dther  story, 
III— 9  through  I 11—15  will  not  require  any  adjustments  which  will  cause 
the  operator  to  get  closer  than  two  feet  to  the  direct  beam.     If  a 
suitable  beam  catcher  is  used,  the  radiation  levels  in  the  occupied 
areas  would  probably  not  exceed  10  mr/hr. 
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4.0  RESULTS 
4.1  Initial  Measurements 

Two  initial  measurements  which  are  made  with  any  neutron  spectrometer 
are  the  rocking  curve  and  the  zero  point  determination.  The  rocking  curve 
is  then  valid  only  so  long  as  the  collimator  openings  and  the  diffracting 
crystal  are  not  changed.  The  zero  point  is  valid  only  so  long  as  the  base 
of  the  instrument  is  not  moved. 

Typical  rocking  curves  taken  with  the  automatic  drive  on  the  crystal 
table  are  shown  in  Fig.  22.  The  time  constant  of  the  ratemeter  was  one 
second  and  the  reactor  power  was  96  kw.  These  curves  were  obtained  at  the 
N.li.  beam  port  at  four  different  diffraction  angles  using  a  (111)  Li F 
crystal.  All  four  were  taken  with  the  crystal  table  rotating  in  the  same 
direction  and  all  yield  a  value  for  the  FWHM  of  the  recking  curve  of  27. 
±  4.  minutes.  Solution  of  equation  17  yields  a  result  of  24  minutes. 
Thus,  it  may  be  assumed  that  the  crystal  and  the  collimator  are  reasonably 
well  matched. 

The  slight  irregularities  in  the  curves  are  caused  primarily  by  slippage 
in  the  drive  belt  from  the  motor.  The  relative  scale  factors  are  given  for 
each  of  the  curves  in  order  for  the  reader  to  obtain  an  estimate  of  the 
magnitudes  of  both  the  diffracted  beam  and  the  background  for  the  different 
diffraction  angles.  Fig.  23  shows  a  typical  zero  point  determination  curve. 
It  was  obtained  at  the  S.E.  beam  port  with  a  reactor  power  of  100  watts. 


These  and  all  following  measurements  were  taken  with  a  h  inch  wide  opening 
on  the  entrance  and  exit  of  the  collimator,  and  a  1  inch  wide  opening  en 
the  front  of  the  detector  shield.  All  openings  were  defined  with  .020  inch 
thick  cadmium  sheet. 
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Fig.  23.  Zero  point  determination  c. 


The  curve  is  shewn  both  with  and  without  the  standard  correction"  for  t!  - 
resolving  time  of  the  counter.  Note  that  the  only  effect  of  this  correc- 
tion is  to  sharpen  the  peak  of  the  curve   and  is  not  really  necessary  in 
.r  to  determine  the  zero  point  location.  The  correction  is  included 
here  only  to  give  an  idea  of  the  order  of  magnitude  of  the  correction  in 
preparation  for  the  determination  of  the  neutron  spectrum  in  section  4.3. 
Another  important  initial  measurement  is  the  reproducibility  of  the 
results.  If  the  results  are  consistent »  suitable  correction  factors  may 
at  least  be  devised  in  order  that  the  true  answers  may  be  obtained.  If 
the  settings  of  the  instrument  are  not  reproducible.,  the  machine  is  nearly 
/orthless. 

The  reproducibility  was  determined  by  selecting  a  set  of  five  differ- 
ent points  in  the  vicinity  of  the  peak  velocity.  Five  determinations  c.r 
each  set  were  made  with  each  point  being  taken  in  cyclic  order.  Since 
the  reactor  power  was  fairly  constant  (constant  resolving  time  correction 
for  each  different  position)  and  the  preset  count  on  the  diffracted  beam 
cons  ant,  only  tae   monitor  counts  are  listed  in  Table  II. 


The  resolving  time  for  the  scaler  and  the  pre  amp  was  experimentally  de- 
mined  by  Ihc   12)   to  be  108  microseconds.  Since  the  BF„  is  operated 
in  the  proportional  region9  its  deadtime  is  much  smaller  thin  the  rest 
of  the  system  and  may  be  neglected.  As  a  results  no  further  investiga- 
tions were  made  along  this  line. 


.'cor  Counts  for  Repeated 
Settings  of  Main  . 


Run.  .  n  ene 


eV 


072  eV         .05,.  eV 


,04C 


„    s>  _ 


2 
5 


98051 

74( 

69527 

.53 

,417 

97C. 

74316 

^31 

72794 

79073 

96472 

;90 

69 :_ 

72195 

7S 

970: 

74415 

u  o 

18 

75325 

69-. 

73". 

79187 

data  taken  at  .0<  the  other 

jrgies.     This  may  be  attributed  to  the  backlash  in  the  and  the 

-king  the  first  setting. 

4.2     Further  Correction  Fa 
4.2.1     Cadmium  Cut-off  Corre: 

Prelir  .ation  of  the  data  d  a  s-c         .  peak  in  the  vicin- 

of  .4  eV.     This  peak  was  finally  attributed  to  the  use  of  c  .  for 

the  definition  of  a  coll       ted  neutron  beam.     In  the  range  from  .1  eV  to  1. 

,  the  total  neutron  cross  section  for  cad:..  .creases  by  nearly 

orders  of  magnitude.     A  correction  factor,  R,  was  then  derived  in  the  fol- 
consider  an  opening  consisting  of  an  unshielde^  area,  A,, 
and  a  cadmium  shielded  area,  A0.     The  ratio,  R,  of  the  neutrons  passing 
through  trc  -.s  is  given  by 
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+  . .,  e 


cc 
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where  z    is  the   macroscopic  total  neutron  cross  section  for  cadmium  and 
t  ^  is  the  thickness  of  the  cadmium  covering  area  A?.  Upon  inspection 
of  the  total  neutron  cross  section  for  cadmium,  it  is  found  that  between 
.25  eV  and  .7  eVs  the  following  is  a  good  approximation: 

3  4 
\   =  (-=f)  .  (26) 

Substitution  of  equation  26  into  equation  25  and  the  introduction  of  the 
proper  numei  ical  quantities  for  the  two  areas  and  the  cadmium  thickness 
yields 

R  = ^— Q2TT  W 

8  +  13  e    E4 

A  graph  of  this  equation  is  shown  in  Fig,  24.  Even  though  equation  25  is 
not  valid  for  energies  less  than  025  eV  or  greater  than  .7  eV,  the  limit- 

j  values  have  nearly  been  reached  at  these  extremes  and  so  the  correc- 
tion factor  may  be  applied  to  all  energies  within  the  range  of  the  instru- 
ment without  any  loss  of  validity  of  the  final  answer. 

The  correction  factor  derived  is  thus  a  correction  for  the  increase 
in  radiating  area.  Although  it  would  seem  that  a  similar  correction  factor 
should  be  applied  to  the  collimator  exit,  no  such  factor  is  needed.  The 
size  of  the  crystal  is  such  that  the  effective  crystal  width  presented  to 
the  collimated  neutron  beam  above  the  cadmium  cut-off  energy  is  less  than 
the  width  of  the  collimator  opening. 

It  will  be  noted  that  in  some  of  the  following  curves,  especially  Fig. 
31,  all  of  the  cadmium  peak  has  not  been  removed.  This  is  attributed  to 
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a  slight  misalignment  of  the  diffracted  neutron  beam.  This  misalic  >d 

..  then  strikes  the  cadmium  slits  of  the  front  of  th  detector  shield. 
Since  the  amount  of  this  misalignment  was  not  known,  no  attempt  was  made 
to  apply  a  suitable  correction. 


.ion  Fa 

e  original  data,    . - 
le  resc .  .  of  the  data  i 

3  correct!  on  )ul 

how \  lied. 

jard  correction  for  resolv"  =   af  a  counter 

+  Rm  Rt>  (- 

:  n     is  the  -ate,  R    is   :  is  the 

o  ^ 

jviously  mentioned,  R.   was  found  tc  :rosecoi 

3  data  collection  and  processing,  equation  28  was 
t  i      ifi 

:    ■  needed  -co  collect  a  givar.  .  counts, 

in  tl  itor   -       nel,  then  t!  -»  Ti  2cJ   -■  z-'-~C:C':  arW  other 

of  counts ,  ...  is 

T.   -  J-x  (29) 

j  reactor  power  is  nearly  constant.     Thus,  only  one  time  need  be 
entered  for  each  di  .     The  .  -  i  resolving  time  correction 

the    .  obtained  was  a^out  1C 

its  frc         rious  Bee 

4.3.1     S.E.  Beamport 
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Figures  25  and  26  show  the  spectra  obtained  for  fuel  loadings  6  and  6' 
using  a  (111)  crystal.  Upon  comparison  of  the  two  spectra ,  it  may  be  seen 
that  the  major  effects  of  changing  between  these  two  core  configurations  is 
to  increase  the  magnitude  of  the  1/E  tail  by  approximately  50%.  Little 
change  is  otherwise  made  in  the  curves. 

Inspection  of  the  original  data  indicated  that  the  countrate  at  the 
peak  of  the  distribution  was  changed  by  4%  with  the  change  in  core  config- 
urations. This  latter  figure  is  not  entirely  meaningful  since  adjustments 
were  made  to  the  crystal  between  these  two  different  runs.  The  nature  of 
these  adjustments  is  such  that  it  is  very   difficult  to  reposition  the  crystal 
so  that  the  diffracted  beam  intensity  is  within  5%  of  that  from  a  previous 
run. 

In  Figures  24  through  36 >  only  a  few  of  the  standard  deviations   are 
shown.  Most  of  those  shown  are  either  in  the  1/E  region  or  close  to  it. 
The  standard  deviations  for  the  remaining  sections  of  the  curves  are  in 
the  neighborhood  of  1-2%  and  are  too  small  to  indicate  on  the  graphs.  Also9 
the  other  uncertainties  for  the  various  correction  factors  are  far  in  excess 
of  those  indicated  in  the  figures.  Besides  the  correction  factors  previously 
mentioned,  it  is  thought  that  the  process  of  off-setting  the  crystal  for 
background  measurements  may  have  contributed  to  the  error  in  the  background 
for  small  diffraction  angles. 


* 

Pertinent  parameters  of  the  experimental  curves  are  given  in  Table  3. 

**These  deviations  are  computed  from  the  measured  values  of  data  channel 
and  monitor  channel  for  both  the  diffracted  beam  and  the  background. 
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For  energies  lower  than  v  ,  the  theoretical  curve  was  obtained  by 
least-squares  fitting  the  normalized  experimental  data  points  to  equation 
6.  The  1/E  tail  is  fitted  usi.'.g  an  eyeball  approach  since  the  data  in 
this  region  are  somewhat  uncertain. 

Interesting  results  were  obtained  for  the  spectrum  (Fig.  27)  taken 
with  the  (200)  crystal.  The  fuel  loading  was  number  6.  This  spectr 
agrees  fairly  well  with  that  shown  in  Fig.  24  down  to  .055  eV.  Below 
that  value,  the  two  spectra  diverge,  with  the  theoretical  curve  lying 
somewhere  betwe'en  the  two  curves.  No  explanation  can  be  given  for  the 
magnitude  of  this  divergence. 


65 


•  ' 


i     r 


I  1  '   l    '     » L 


o 


rz 

o 


UJ 


b 

>s 

a 

C~ 

Q 

<. 

CO 

i 

O) 

c 

( — 

i 

1_J 

c 

c 

-i-> 

• 

o 

Z^ 

a 


CM 


CD 


O 

o 


C  4-> 

o  c 

S-  CD  -— " - 

+■>  S-  LU 

O  3    ^ 


.  Beam  Po. 

Figure  30  show  the  spectra  obtained  usi 

5" ,  re  ht  be  e       :    sd,  lit...  -  di  Ffer 

.3d  between  the  first  two  fuel  loadings.  Howevei 

magnitude  of  the  1/E  tail  by  nearly  a  factor  of 

:rates  obtain:     fuel 
However.,  t    .  ..  .  countr       uel 
lc^  :er  than  t        i  -rs.  Here, 

made  between  the  tnree  runs. 
Certain  pertu  :inity  of  .01  eV  are  much  more  pron  lent 

in  these  spectra  than  for  the  spectra  from  the  S.E     .  As  will  be  shown  in 

;ion  -. .     ;e  perturbations  are  causec     2  presence  of  a  size. 

of  all     in  the  neutron  path. 
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4.3.3 Beamport 

Figures  31,  32  and  33  show  the  spectra  obtained  using  the  standard 
fuel  loading.  This  series  of  measurements  was  conduci-v.  in  order  to  de- 
li ne  the  effectiveness  of  the  alignment  procedures.  After  the  spectrum 
shown  in  Fig.  31  was  obtained,  the  entire  instrument  was  moved  to  the  N.E. 
port  and  several  spectra  taken.  The  spectrometer  was  then  returned  to  the 
N.W.  port  and  realigned.  The  spectra  shown  in  Figures  32  and  33  were  then 
obtained. 

Inspection  of  the  original  data  indicates  that  the  peak  countrates 
were  within  about  10%  of  each  other  for  the  three  different  runs.  The  data 
in  Fig.  32  are  somewhat  suspect.  Shortly  after  these  data  were  taken,  the 
scaler  on  the  monitor  channel  was  found  to  be  oscillating  intermittently. 

Fig.  34  shows  the  spectrum  obtained  when  a  piece  of  aluminum  2  1/3 
inches  thick  was  placed  in  the  collimator  shielding  hole.  The  peaks  (24) 
in  the  cross  section  around  .01  and  .015  eV  are  especially  pronounced.  The 
effect  of  this  thickness  of  aluminum  was  to  reduce  the  maximum  countrate  by 
a  factor  of  two. 

During  the  realignment  procedure,  it  was  noticed  that  the  longitudinal 
axis  of  the  collimator  was  not  horizontal.  The  entire  collimator  was  rotated 
in  the  beam  tube  180  degrees  to  ensure  that  this  effect  was  not  caused  by  a 
warped  collimator.  Since  this  procedure  did  not  correct  the  misalignment, 
it  was  assumed  that  the  beam  tube  was  at  fault.  No  similar  difficulties 
e  exDerienced  at  either  the  S.E.  or  the  N.E.  port. 
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Table  III.  Pertinent  Values  from  Theoretical  Curves 
Fitted  to  the  Beamport  Leakage  Spectra 


Beam- 
port 

Date    Crystal 

Fuel 
loading 

Kl3 

rsf 2  x  io-9) 

cm 

D 
V 

mp 
(m/sec  x  103) 

vc 
(m/sec  x  103) 

S.E. 

6/16/66 

(111) 

6 

.1416 

3.098 

b.o4 

6/17/66 

(200) 

6 

.1007 

3.342 

6/29/S6 

(111) 

6 ' 

.1416 

3.161 

0.11 

7/5/66 

(111) 

6 

.1332 

0  .  1  uj 

5.88 

7/6/66 

(111) 

6„. 

.1332 

3.177 

6.11 

7/6/66 

(111) 

6" 

.  1 308 

3.228 

5.79 

7/1/66 

(111) 

6 

.2125 

2.774 

5.53 

7/13/66 

(HI) 

6 

.1928 

2.841 

5.71 

7/14/66 

(111) 

6 

.1988 

2.971 

5.71 

a,, 


iaximum  value  in  corrected  experimental  data  has  been  normalized  to  unity. 
}To  obtain  v  ,  multiply  these  values  by  /2/3. 


'No  interception  point  could  be  accurately  determined. 


4.4  Measurement  of  the  Energy  Dependence  of  the 
Total  Neutron  Cross  Section  for  Beryllium 

.1  neutron  cross  section  for  beryllium      is  shown  in 
.  35.  There  is  fair  agreement  for  the  two  different  runs  with  the 
published  values  (25)  above  .01  eV.  The  peak  at  .02  eV  is  even  better 
defined  than  the  published  value.  Below  .  01  eV,  the  agreement  is  not 
so  goo  .  This  may  be  readily  attributed  to  the  predominance  of  the 
third  order  diffraction  pattern  in  this  region. 

Fig.  35  shows  the  variation  of  count  rate  versus  varying  thickness 
of  beryllium  for  different  neutron  energies.  The  nonlinearity  of  these 
curves  also  points  out  the  effect  of  the  higher  order  diffraction  patterns 
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5.0  CONCLUSIONS 

The  measurements  contained  within  this  work  have  shown  that: 

1.  The  KSU  neutron  spectrometer  may  be  operated  safely,  provided 
it  is  operated  according  to  the  instructions  contained  in 
section  3.4.  If  these  instructions  are  faithfully  followed, 
the  maximum  permissable  exposures  contained  in  10CFR20  will 
not  be  exceeded. 

2.  Core  configurations  may  be  chosen  for  beamport  experiments  in 
such  a  way  so  as  to  yield  the  maximum  possible  signal  to  noise 
ratio. 

3.  The  total  neutron  cross  section  for  Be  has  been  measured  and 
found  to  be  in  reasonable  agreement  with  published  values. 
Measurement  of  this  cross  section  has  shown  the  KSU  neutron 
spectrometer  to  possess  wery   good  resolution. 

4.  An  instrument  has  been  constructed  which  will  allow  undergrad- 
uate students  at  the  senior  level  to  measure  the  thermal  neutron 
spectrum  from  the  KSUTMII  and  obtain  reasonable  results  within 
an  afternoon.  Provided  the  instrument  has  been  previously 
aligned,  these  results  will  be  reasonably  accurate. 

5.  There  is  a  definite  perturbation  of  the  thermal  neutron  spectrum 
caused  by  the  presence  of  a  considerable  thickness  of  aluminum 
near  the  entrance  of  the  beam  tubes.  At  the  N.E.  beam  port,  a 
certain  perturbation  may  also  be  noticed  from  the  graphite  dummy 
elements  for  certain  core  configurations. 
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6.  The,  slits  which  define  the  collimated  neutron  beam  should  be 
constructed  of  material  other  than  cadmium  in  order  to  better 
col  lima te  the  neutron  beam  in  the  epi cadmium  region  and  to 
reduce  the  background. 

7.  Improvements  should  be  made  on  the  crystal  table  assembly, 
especially  in  that  portion  which  provides  for  rotation  of  the 
crystal  about  a  vertical  axis. 

8.  Measurement  of  the  total  neutron  cross  section  for  beryllium 
indicates  that  the  energy  calibration  of  the  instrument  is  in 
error  by  7%   at  .021  eV.  Further  studies  should  be  made  to 
determine  the  calibration  error  as  a  function  of  position. 
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6.0  SUGGESTIONS  FOR  FURTHER  WORK 

There  are  a  number  of  improvements  which  may  be  made  in  the  present 
system.  One  of  the  most  important  of  these  is  the  design  and  construction 
of  better  collimating  slits.  A  rough  measurement  showed  that  nearly  90% 
of  the  background  consisted  of  epicadmium  neutrons.  A  better  slit  system 
could  reduce  the  background  by  at  least  a  factor  of  two. 

If  much  work  is  to  be  done  in  the  lower  energy  region,  a  Soller  slit 
system  inserted  into  the  collimator  could  more  than  double  the  present 
countrate  without  reducing  the  resolution. 

The  start  of  a  second  axis  has  been  provided.  Because  of  lack  of  time, 
it  was  not  completed.  The  experience  obtained  with  the  main  support  ring 
indicated  that  the  secondary  support  ring  may  have  to  be  removed  and  another 
means  of  support  found.  This  procedure  would  also  remove  a  great  deal  of 
weight  as  well  as  gain  some  extra  height  for  the  second  gearbox. 

Another  improvement  would  be  the  substitution  of  a  very  thin  ion  chamber 
placed  directly  in  the  neutron  beam  in  place  of  the  present  monitor  probe. 
This  would  guard  against  flux  tilts  within  the  core. 

The  above  improvements  would  allow  more  work  to  be  done  on  the  higher 
order  diffraction  patterns  and  refinement  of  the  Be  cross  section  measurements 
Also,  they  would  greatly  increase  the  versatility  of  the  instrument. 

Another  interesting  experiment  which  is  suggested  is  the  study  of  the 
thermal  neutron  spectrum  as  a  function  of  time  during  a  reactor  pulse.  Ob- 
viously, this  would  have  to  await  a  pulsing  license.  It  is  suggested  that 
if  this  experiment  is  performed,  a  stepped  hole  be  cut  in  the  wall  of  the 
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neutron  generator  room  directly  opposite  the  fast  beam  port  in  order  to  use 
the  entire  room  as  a  beam  catcher.  Naturally,  this  hole  could  be  plugged 
during  operation  of  the  neutron  generator. 
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9.0  APPENDICES 

9.1  Appendix  A:  Beam  Tube  Measurements 

9.1.1  Physical  Dimensions  of  Beam  Tubes 

The  physical  dimensions  of  the  beam  tubes  which  were  required  for  the 
collimator  design  are  given  in  Tables  IV  and  V.  These  dimensions  may  be 
used  in  the  design  of  other  collimators  or  for  determining  the  clearance 
between  the  end  of  the  present  collimator  and  the  end  of  the  beam  tube  for 


use  in  future  experiments. 


Reactor  Shield 
h^r. A 


Reactor  Tank  Wall 


c 


Shadow  Sh 


Fig.  37.  Cross  section  of  typical  beam  tube 


Table  IV.  Longitudinal  Dimensions  of  Beam  Tubes 


Beam  tube 


Dimension  (in  inches) 
ABC 


S.E. 

N.E. 
N.W. 

s.w. 


64 

30-1/4 

112-1/8 

64-1/4 

29-7/8 

110-15/16 

64-1/4 

30-1/4 

m-i/4 

64-1/4 

29-7/8 

110-15/16 

37 


General  Atomic  (19)  gives  the  inside  diameters  for  all  four  beam 
tubes  as  follows. 

Table  V.  Inside  Diameters  of  Beam  Tubes 


Dimension  Value  in  inches 

D  8  (nominal) 

E  6.065 

F  6.315 

9.1.2  Summary  of  Radiation  Survey  Measurements  from  Open  Beam  Ports 

9.1.2.1  Direct  Beam 

These  measurements  are  included  to  allow  the  reader  to  gain  an 
idea  of  the  radiation  hazard  involved  with  the  direct  beam  from  a  reac- 
tor beam  port.  Measurements  are  given  for  a  completely  open  beam  port 
as  well  as  a  collimated  beam. 

Dose  rate  measurements  taken  by  placing  dosimeters  in  the  colli- 
mated radiation  beam  yield  results  which  are  in  reasonable  agreement 
with  those  in  Table  VI.  For  this  series  of  measurements,  two  neutron 
and  gamma  sensitive  dosimeters  and  two  gamma  sensitive  dosimeters  were 
held  with  tongs  in  the  collimated  radiation  beam.  An  attempt  was  made 
to  time  the  irradiation  so  that  the  dosimeters  would  reach  mid-scale 
regardless  of  reactor  power  level.  Two  additional  neutron  and  gamma 
sensitive  dosimeters  were  irradiated  at  the  maximum  power  level.  These 
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Table  VI.  Measurements  of  Direct  Radiation  Beams 
from  Completely  Open  Beam  Ports 


Beam 

Power 

Neutrons 

Gamma  (mr/hr) 

tube 

level 
(watts) 

Hurst  detector      Nemo 
Fast*     Slow** 
(cps)     (cps) 

(mrem/hr) 

S.E. 

1 

5.6 

3 

5 

25 

10 

15 

77 

30 

50 

280 

100 

150 

860 

300 

500 

— 

800 

N.E. 

1 

13       400 

75 

24 

5 

160      1000 

380 

70 

10 

780 

100 

N.W. 

1 

11 

10 

5 

52 

30 

15 

130 

70 

♦Multiply  by  1  to  obtain  mrem/hr. 
**Multiply  by  .08  to  obtain  mrem/hr. 
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measurements  are  shown  in  Fig.  38.  It  would  appear  that  the  two  gamma 
sensitive  dosimeters  were  not  irradiated  for  the  proper  length  of  time 
at  96  kw. 

In  order  to  compare  the  dosimeter  measurements  with  those  in  Table 
VI,  a  correction  must  be  made  for  the  radiating  area  as  seen  by  the  de- 
tector. The  magnitude  of  this  correction  has  been  calculated  to  be  40 
(26).  If  the  total  dose  rate  from  the  three  beamports  as  given  in  Table 

VI  is  extrapolated  to  a  reactor  power  level  of  100  kw  and  then  divided 
by  the  radiating  area  correction  factor,  the  dose  rates  shown  in  Table 

VII  are  obtained. 

Table  VII.  Collimated  Beam  Total  Dose  Rates  for 
a  Reactor  Power  Level  of  100  kw 


Beamport  Dose  rate  (rem/hr) 

S.E.  19 

N.E.  220 

N.W.  30 

If  calculations  are  made  of  the  separate  contributions  to  the  dose 
rate  from  both  the  gammas  and  the  neutrons,  it  is  seen  that  while  the  gamma 
dose  rate  agrees  to  within  50%  of  the  dosimeter  data,  the  corresponding 
value  for  the  neutron  dose  rate  does  not.  This  is  rather  easily  explained 
since  the  dosimeters  which  were  used  are  only  sensitive  to  thermal  neutrons 
while  the  NEMO  measures  dose  rates  for  all  energy  neutrons. 
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9.1.2.2  Scattered  Radiation  Dose  Rates 


Positive  control  procedures  make  the  possibility  of  operating  the 
KSUTMII  reactor  with  accidentally  opened  beam  tubes  vani shingly  small. 
In  order  to  gain  an  idea  of  the  order  of  magnitude  of  the  radiation 
levels  involved  in  such  an  occurence,  it  v/as  deemed  desirable  to  operate 
the  reactor  in  such  a  condition.  The  shielding  plugs  were  completely 
removed  from  the  N.E.  beam  tube  and  the  reactor  bay  evacuated.  The 
power  level  was  raised  in  steps  as  determined  by  the  radiation  safety 
officer  and  the  reactor  supervisor.  Measurements  were  then  carefully 
taken  at  selected  points  within  the  reactor  bay.  These  measurements  are 
given  in  Table  VIII. 

Because  of  the  limited  space  between  the  N.E.  beam  port  and  the 
shielding  wall  for  the  neutron  generator  room,  it  was  decided  to  investi- 
gate the  possibility  of  using  this  wall  for  a  beam  catcher.  Since  a  fast 
neutron  scattering  experiment  was  already  in  progress  at  this  beam  port, 
the  graduate  student  involved  was  prevailed  upon  to  calculate  a  dose  rate 
from  his  data.  The  particular  geometry  for  his  experiment  is  given  in  Fig. 
39. 


Concrete  slab 


Detector 
Fig.  39.  Experimental  arrangement  for  the  determination  of  the  dose  rate  from 

fast  neutrons  backscattered  from  a  concrete  slab 
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Table  VIII.  Radiation  Levels  in  Reactor  Bay  with 
Unplugged  Collimator  in  N.E.  Beamport 


Position 

Power 

(kw) 

.1 

1. 

10. 

25. 

Gamma  dose  ina 
direct  beam 

2. 

10. 

40. 

2  feet  north 
of  beam 

1.5b 

64. 
100. 

600. 

- 

neutron  gener- 
ator console 

3.5 
.75 

50. 
10. 

400. 
200. 

- 

east  corner  of 
loading  dock 

.46 
.1 

3.8 
.12 

50. 
.9 

100. 
16. 

desks  in  N.W. 
corner 

.7 

.025 

5.5 
.08 

30. 
1.9 

94. 
13. 

control 
room  wall 

.17 
.08 

.55 
.06 

25. 
.9 

40. 
5.5 

22  foot  level 
N.E.  corner 

.32 
.1 

3. 
.13 

180. 
1.15 

140. 
65. 

Secondary 
water  tank 

1. 
.025 

3. 

.05 

18. 
3.5 

52. 

8.5 

Barnstead 
still 

.6 
.025 

.5 
.025 

1.8 
.13 

4.8 
.9 

N.E.  corner 
of  bay,  at 
wall 

- 

1.2 

no. 

500. 

These  readings  only  in  R/hr. 

}Top  reading  is  neutron  dose  rate  in  mrem/hr. 
Bottom  reading  is  gamma  dose  rate  in  mr/hr. 
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The  results  thus  obtained  (27)  are  given  in  Table  IX.  One  of  these  posi- 
tions was  checked  with  a  film  badge.  The  low  dose  obtained  yielded  incon- 
clusive results.  As  a  result  of  the  values  obtained  in  Table  IX,  it  was 
decided  to  use  a  beam  catcher.  The  beam  catcher  finally  used  for  this 
beam  port  consisted  of  the  detector  shield  used  to  bbtain  the  results  in 
Fig.  39. 

Table  IX.  Fast  Neutron  Dose  Rate  Measurements  of  Radiation  Scat- 
tered from  a  Concrete  Slab  (Reactor  Power  is  10  kw) 

a  Fast  neutron  dose  rate 

(mrem/hr) 

19  21.6 

33  20.4 

44  15.4 

52  13.4 

60  11.8 

67  10.4 
73  8.28 

80  4.60 

Table  X  shows  a  survey  taken  in  the  vicinity  of  the  neutron  spectrom- 
eter at  the  S.E.  beam  port.  Surveys  have  been  carried  out  for  the  operation 
at  N.E.  and  N.W.  beam  ports  as  well.  Since  they  all  yield  nearly  the  same 
results,  only  the  one  is  included  here. 
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Table  X.  Survey  around  Neutron  Spectrometer 
at  S.E.  Beamport 


Power 

Gamma 

Fast 

Slow 

Location 

level 

(mr/hr) 

neutrons 
(cps) 

neutrons 
(cps) 

10  kw 

.5-4 

at  crystal  and  beam  exit 

15 
2 

20 
<1 

600 
<1 

front  of  beam  catcher  (20* 

from  port) 
back  of  beam  catcher 

<1 

<1 

<1 

periphery,  at  rope 

50  kw 

3-5 

2 

20 

at  face  of  beamport 

3-5 

2 

40 

Rowland's  position3 

70 

30 

1000 

front  of  beam  catcher 

7 

<1 

<1 

back  of  beam  catcher 

100  kw 

no 

40 

off  scale 

front  of  beam  catcher 

13 

2 

7 

back  of  beam  catcher 

5 

3 

40 

Rowland's  position3 

2 

- 

10 

periphery 

1 

2 

10 

at  instrumentation  table 

aThis  position  is  that  occupied  by  the  experimenter  when  making  adjustments 
to  the  crystal  as  per  section  3.4. 

Nuclear  Chicago  Ion  Chamber  was  used  for  gamma  detection.  The  n,  y 
reaction  may  have  helped  to  increase  the  gamma  dose  rate  level. 

The  RCL  counter  (N.E.  No.  585)  was  used  for  neutron  detection.  To 
convert  fast  neutrons  to  mrem/hr,  multiply  by  one  (1);  to  convert  slow 
neutrons  to  mrem/hr,  multiply  by  .08. 

John  P.  Lambert 
Radiation  Safety  Officer 
12/1/65 
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9.1.3  Cadmium  Ratio 

Very  simple  measurements  of  the  cadmium  ratio  were  made  at  the  S.E. 
and  N.E.  beamports.  The  values  obtained  were  47  and  17,  respectively. 
These  measurements  were  taken  only  as  a  matter  of  curiousity,  but  should 
be  accurate  to  within  30%. 
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9.2    Appendix  B:  Activation  of  Surplus  Lead  Samples 

Two  lead  samples  were  obtained  through  State  Surplus  in  Topeka.  A 
one  minute  irradiation  was  found  to  be  sufficient  for  the  .6  gram  samples 

used.  The  reactor  power  level  was  100  kw  and  the  flux  at  the  point  of 

11     2 

irradiation  (rotary  specimen  rack)  was  approximately  7  x  10   n/cm  -sec. 

The  results  are  given  in  Table  XI. 

Two  slightly  larger  samples  (obtained  from  the  same  sources)  were 
irradiated  in  the  AECL  GAMMA  CELL  for  nearly  24  hours.  No  evidence  of 
photoactivation  was  observed. 

A  saturation  activity  test  was  run  on  a  .25  gram  sample  of  the  surplus 
lead.  The  power  level  was  100  kw  and  the  irradiation  took  place  first  in 
the  rabbit  and  then  in  the  rotary  specimen  rack.  Saturation  seemed  to 
occur  after  five  minutes.  This  test  was  only  continued  for  40  minutes. 
The  saturation  activity  without  the  beta  shield  was  75  mr/hr.  This  reading 
was  taken  within  five  minutes  after  the  sample  was  removed  from  the  reactor. 

From  the  results  shown  in  Table  XI,  it  may  be  seen  that  there  are  at 
least  two  different  isotopes  involved.  One  of  these  is  short  lived  (5-10 
minute  half  life  or  less)  and  is  probably  a  beta  or  a  weak  gamma  emitter. 
One-eighth  inch  of  aluminum  plus  1-1/2  inches  of  air  were  sufficient  to 
reduce  the  saturation  activation  reading  of  75  mr/hr  to  4  mr/hr. 

The  other  isotope  apparently  has  a  half  life  of  about  2  1/2  days  and 
is  also  a  beta  or  very  weak  gamma  emitter. 

From  the  half  life  data,  it  is  obvious  that  the  saturation  activity  test 
was  not  continued  long  enough  to  produce  a  complete  result.  It  does  tend  to 
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indicate,  however,  that  the  half  life  of  the  short  lived  isotope  is  more  on 
the  order  of  1-2  minutes  than  5-10  minutes.     Consultation  of  appropriate 
half  life  data  indicates  that  the  impurity  is  likely  to  be  antimony. 
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Table  XI.  Decay  History  of  Activated  Surplus  Lead  Samples 


K  K 

Date     Time     Sample     Readings  outside    With  beta    Without  beta 

rabbit  without      shield       shield 
beta  shield 


8/2/65 

2:10pm 

surplus 

50a 

— 

—  — 

2:12pm 

site 

>200 

— 

~ 

3:15pm 

surplus 

11.5 

6.5 

130 

3:15pm 

site 

11.5 

9 

120 

7:15pm 

surplus 

10.5 

6 

125 

7:15pm 

site 

9 

6 

125 

8/3/65 

2:55pm 

surplus 

5.5 

5 

no 

2:55pm 

site 

6.5 

5.5 

100 

8/4/65 

3:50pm 

surplus 

5.5 

4.5 

105 

3:50pm 

site 

4.5 

4 

90 

8/15/65 

7:40pm 

surplus 

.5 

.6 

5.5 

7:40pm 

site 

.7 

.6 

5 

aAll  readings  in  mr/hr  taken  with  a  standard  geiger  survey  meter. 
Sample  was  placed  next  to  the  detector. 
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9.3  Appendix  C:  Collimator  Effectiveness  and  Safety 

The  measurements  contained  within  this  appendix  show  that  the  collimator 
provides  excellent  collimation  for  thermal  neutrons.  They  also  prove  that 
the  collimator  shielding  plugs  allow  the  reactor  to  be  operated  at  full  licen- 
sed power  with  a  plugged  collimator  in  place  without  producing  any  significant 
radiation  hazard.  It  has  also  been  demonstrated  that  the  absence  of  shielding 
plugs  in  the  collimator  may  readily  be  detected  by  the  reactor  operator  even 
though  the  outer  beamport  door  is  closed <,  provided  the  RAM  probes  for  the 
beamports  are  placed  in  the  proper  positions. 

The  effectiveness  of  the  collimator  for  thermal  neutrons  was  determined 
by  moving  a  small  vertical  BF^  probe  horizontally  across  the  collimated  neut- 
ron beam  at  a  constant  velocity.  The  data  were  accumulated  by  a  TMC  GAMMA 
SCOPE  operating  in  the  multiscaler  mode.  The  timing  of  the  multiscaler  was 
one  channel  per  second  and  the  traversing  mechanism  traveled  at  one  inch  per 
minute.  Although  the  diameter  of  the  BF~  probe  was  .3  inch9  no  correction 
for  this  was  made  in  the  final  data  since  only  an  estimation  of  the  flux 
profile  was  desired.  Fig.  40  gives  results  for  probe  to  snout  distances  of 
two  and  six  feet.  These  results  agree  quite  well  with  values  calculated 
using  geometrical  considerations. 

Table  XII  shows  radiation  measurements  indicating  the  effectiveness  of 
the  collimator  shielding  plugs.  The  collimator  was  placed  in  the  N.E.  beam- 
port  since  this  beamport  provides  the  maximum  gamma  and  fast  neutron  dcs^ 
rates.  With  both  collimator  plugs  in  place  and  the  safety  shutter  and  the 
outer  door  closed,  the  maximum  dose  rate  around  the  beamport  -is  less  than  2 
mr/hr. 
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Table  XIII  shows  the  RAM  readings  as  a  function  of  reactor  power 
level.  To  obtain  these  measurements,  the  collimator  was  placed  in  the  S.E. 
beamport  without  its  internal  shielding  plugs  and  the  outer  door  of  the 
beamport  locked.  The  RAM  probe  was  hung  on  the  outside  of  the  door  and 
the  readings  in  the  control  room  noted  as  the  reactor  was  brought  up  to 
power.  At  the  maximum  power  in  each  series,  a  rough  survey  was  taken  of  the 
immediate  area  near  the  door.  Two  different  series  were  taken,  one  with  the 
safety  shutter  open  and  one  with  it  closed.  Also,  the  effect  of  placing 
the  RAM  probe  in  different  positions  on  the  door  was  noted.  As  a  result 
of  these  measurements,  it  is  suggested  that  all  beamport  RAM  probes  be  hung 
on  the  beamport  doors  as  indicated  by  position  1  in  Fig.  41 .  This  position 
placed  the  detector  closest  to  the  center  of  the  beam  tube. 


Fig.  4(1.  Different  positions  for  placing  the 
RAM  probe  on  beamport  doors 
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Table  XII.  Radiation  Measurements  of  Plugged  Collimator 


Condition  of  safety             NEMO  Geiger  survey  meter 

shutter  and  beamport  door        neutron  gamma  dose  rate 

dose  rate  mr/hr 
mrem/hr 


both  open  4.0  8.5 

both  closed  .75  .15 


Table  XIII.  RAM  Readings  for  Unplugged 
Collimator  in  S.E.  Beamport 


RAM  probe 

Power 

RAM 

reading  with 

position3 

level 

safety 

shutter 

(kw) 

and 

door  closed 

1 

1 

10 

50 

90 

.1 

1 
10 
40 

1 

20 

125 

250 

2 

40 

3 

40 

RAM  reading  with   Area 
only  door  closed  survey 

results0 


8 
1.2 
15 
220 

920        >100 
500 
300 


aRefer  to  Fig.  41. 


DUnless  otherwise  stated,  all  readings  are  in  mr/hr, 
cCombine  neutron  and  gamma  dose  rates  in  mrem/hr. 
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9.4    Appendix  D:     Determination  of  the  Mosaic  Spread 
of  the  Dif     acting  Crystals 

The  mosaic  spread  of  the  diffracting  crystals  was  determined  by  the 
conventional  x-ray  diffraction  method.     A  calcite  crystal  was  used  to 
select  the  cobalt  K     line  from  a  suitable  x-ray  tube.     This  monobromated 
beam  was  then  directed  upon  the  particular  Li F  crystal  under  investigation. 
A  geiger  counter  was  placed  at  the  appropriate  position  and  the  LiF  crystal 
was  then  slowly  rotated  through  the  Bragg  angle  at  a  constant  angular 
velocity.     The  results  obtained  are  summarized  in  Table  XIV. 

As  one  may  see  from  Fig.   4  2,  it  is  not  surprising  that  there  is  sue;, 
wide  disagreement  between  the  different  values  in  Table  XIV.     The  irregu- 
larities in  the  curve  arise  as  a  result  of  the  small   area  actually  being 
measured.     It  would  have  been  more  proper  to  use  a  neutron  beam  of  large 
enough  dimensions  so  that  the  entire  face  of  the  crystal  would  be  illuminated, 
However 9  the  equipment  was  not  available  to  do  this  to  the  required  accuracy. 
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Table  XIV.  Measured  Mosaic  Spreads  and  Peak  Inten- 
sities of  Li F  Diffracting  Crystals 


Crystal 

Surface 

n(FWHM) 

Relative 

(minutes) 

Peak  Intensity 
(cpm) 

(100)  (#2) 

cleaved 

3-1/4 

.78  x  200 

(100)  (#2) 

ground 

8-3/8 

.83  x  500 

(100)  (#2) 

ground 

9-3/4 

.70  x  500 

(100)  (#4) 

cleaved 

1-1/2 

.67  x  500 

(111)  (#3) 

sawed 

12-3/4 

.91  x  500 

(11.1)  (#D 

sawed 

16-1/4 

.76  x  500 

(HI)  (#1) 

sawed 

12-1/2 

.62  x  500 

(111)  (#D 

sawed 

17-1/2 

.80  x  500 
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9.5  Appendix  E:  Explanation  of  Computer  Codes 
(J.5.1   REFLECTIVITY 


* 


The  REFLECTIVITY  computer  code  was  written  to  numerically  integrate 
equation  15  as  a  function  of  diffraction  angle  and  mosaic  spread.  Since 
this  is  a  one-shot  program,  no  input  data  is  necessary.  Instead,  the 
values  for  the  mosaic  spread  and  the  diffraction  angle  are  contained  with- 
in the  program.  The  diffraction  angle  is  determined  in  a  DO  loop  which 
determines  12  different  angles  starting  with  five  degrees. 

The  computer  code  consists  of  three  parts.  The  first  part,  REFLECT, 
performs  the  following  tasks: 

1.  computes  the  values  of  constants  involving  the  mosaic  spread 
and  the  diffraction  angle, 

2.  computes  the  appropriate  limits  of  integration, 

3.  sets  the  necessary  internal  switches  for  the  subprograms,  and 

4.  prints  the  final  answers.    Each  line  of  answers  contains  the 
mosaic  spread,  diffraction  angle,  and  the  associated  integrated 
reflectivity. 

Table  XV  describes  the  important  variables  contained  in  REFLECT.  As  shown, 
REFLECT  contains  the  constants  for  a  Li F  (200)  crystal.  For  a  Li F  (111) 
crystal,  the  following  change  in  statement  1  +  13  should  be  made:  V  = 
.00604*  (the  rest  is  the  same).  Also,  ETAA  should  be  changed  to  reflect 
the  different  value  of  the  mosaic  spread. 


*A11  computer  codes  in  this  work  are  written  in  FORTRAN  IV  for  the  PR-155 
interpretive  system  on  the  IBM  1410  computer. 

For  the  proper  input  and  output  formats  for  all  programs,  refer  to  the 
appropriate  program  listing. 
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Table  XV.  Explanation  of  Selected  Variables  in  REFLECT 

Symbol  Explanation 

A  A  as  defined  in  section  2.5.1 

ANSWB  R  as  defined  in  equation  15 

ETAA  mosaic  spread  in  seconds 

GA  lower  limit  of  integration 

GB  upper  limit  of  integration 

THETAA  diffraction  angle  in  degrees 

V  V  as  defined3  in  section  2.5.1 

The  values  for  the  structure  factor  used  in  calculating 
this  variable  are  experimental  values  (28)  rather  than 
the  theoretical  values. 

Table  XVI.  Explanation  of  Selected  Variables  in  DUMMY 

Symbol  Explanation 

A  A  as  defined  in  section  2.5.1 

T  T  as  defined  in  section  2.5.1 

V  V  as  defined  in  section  2.5.1 


. 


The  function  subprogram,  DUMMY,  calculates  the  value  of  equation  15 
at  the  points  determined  by  the  numerical  integration  subroutine.  The  im- 
portant variables  of  DUMMY  are  explained  in  Table  XVI. 

GRAL  is  a  numerical  integration  subroutine  obtained  from  the  KSU 
computing  center.  A  complete  explanation  of  the  subroutine  is  available 
from  this  source. 

As  shown,  the  program  requires  about  40  minutes  to  run.  The  results 
obtained  for  the  two  different  crystal  planes  are  given  in  Appendix  F. 
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MON$,$      JOB   REFLfCTlVlTY   JIM  ROWLAND 

MON$$      CCMT  15*5PAGES*R0WLAND»NUCLEAR  ENGG   

MCN$$      ASGM  MJB,12 

VON'**      ASGN  MGC»16  "' "~ '  ~ "  " 

MON$$      MODE  GO»TEST 

MON$$      EXEO  FCRTRAN»»#18»»f »REFLECT 
EXTERNAL  DUMMY 

COMMON  A»V  ■  •  '  - -r- 

1  FORMAT (3E 18. 8) 

GA  =  1. 

GA^  =  0.  .  - 

FTAA=9. 

FTA=FTAA*?. 141 5926/1 0800. 

THETAA=0. 

GRE=.OO01 

THETAA=THETAA+5.  ~   

THETA=THETAA*3. 1415926/180.  . 

A=2.69/.3867  

F=SIN(THETA) 
V=.00256*.38  67*F*F/(ETA*SIN(2.*THETA) ) "*" 

GR=1.+?.*V 

AM.S=GRAL(DUMMY»GA»GR»GAF»GRF»IGUA) 

ANSB=FTA*ANS*2.828 
10  wpiTF(3,l  )PTAA,THFTAA»ANSS  "" ' 

CALL  EXIT 

END  

MON$$      EXEO  FORTRAN* »»18 

FUNCTION  DUMMY(T) 

COMMON  A*V  ■    •  .. 

HT=SORT(T) 

T0P=EXP(2.*A*HT)-1. 
C*$*$$  DENOM  OF  LOG  CHANGED  TO  AVOID  TAKING  LOG(U/0)  ■  ; 

B0T=(T+1. )*TOP*ALO<  (2.*V/( T-. 999999) ) +4.*HT* ( TOP+2 • ) 
"'  DUMMY  =  TOP/BOT 

^FTURN 

END 
.   . ,MON$$      FXEO  F0RTRAN»»*18 

FUNCTION  GRAL(DUMM> »GA *GB »GAE *GRE ♦ I GUA ) 
C      THIS  FUNCTION  COMPUTES  THE  INTEGRAL  OF  THE  FUNCTION  GIVEN 
C     IAS  A  FUNCTION  TYPE  SUBPROGRAM  IN  THE  CALLING  STATEMENT 
C     2FR0M  GA  TO  GB  WITH  AN  ABSOLUTE  ACCURACY  OF  GAE  AND  A 
C     3RELATIVE  ACCURACY  OF  GRE. 

C      IGUA=0  FOR  A  NORMAL  I NTERGRAT ION*  IGUA=1  FOR  AN 
C     1FXTENSI0N  OF  A  PREVIOUS  INTFPGRAT ION. " 

DIMENSION  GW(5) 

1041  FORMAT! IX, 44HB  IS  NOT  GREATER  THAN  PRESENT  VALUE  CF  GXMOM ) 

1042  F0RMAT(1X*?3HB  IS  NOT  GREATER  THAN  A) 

1043  F0RMAT(1X,14HH  IS  TOO  SMALL *5X »6HXMCM  = » 1PE18 .7»5X» 3 
?HH  'e.FIA.71  - 


no 


-  C  PFGIN  •   **• — 

•GPAL    =    0.0 

1000  IF( IGUA) 1( 01 .     1001*     1003  -     -- 

1001  IF(GB-GXMCM) 1002*     1002*    1009 
10  02    WRI TE(3»1041 ) 

1C44    RFTURN 

100?  IF(CR-GA) 1004*  1004*  1005  •— -"  " 

1004  WRI TF<3»1042) 
r,r  TO  104  6 

1005  GOAL  =  0,0 
GXMCM  =  GA 
GMAMI  =  0.0 

—  GXMAX  =  GA  '       ~  " r — ~* 

GHXMA  =0.0  • 

155 1 GM  =  0.0 

GFXMO  =  DUMMY (GA) 

GHSUG  =  0.25*(GC-GA)   ' 
1009  IGSWA  =  0 

~GHMMI=u.025*(GB-GXMCM)*GRE  

C      STOP 
1012  IF(GH<;i!G-GHMMI  )1013»  1014,  1014 
101?  GHSUG  =  GH^MI 
1014  GH  =  GHSUG 

IGSWP  =  0 


IF(GXMCM+4.*GH-GB)1016»  1015V  1015      "  ". 

1015  GH  =  0.25*(GR-GXMCM) 
TGSWA  =  1 

C      YFT  CME 

1016  GX  =  GXMCM  "-    ~~ " 
GW( 1 )  =  GFXMC 

00  1020   TGK=2,5  "  " TTm " 

IF( IGK-5H019,  1017,  1017 

1017  TF( TG^WA) 1019,  1019,  1018 

101 8  gw( IGK)  =  DUMMY ( GE ) 
GC  TO  K  20 

1019  GX  =  GX+GH 

GW(IGK)  =  DUMMY(GX) 

1020  CONTINUE 

GMID  =  ( 14.*(GW< 1 )+GW(5) )+64.*(GW< 2)+GW<4) )+24«*GW<3) )/45r 
GDELT=ABS( (2.*(GW( 1 )+GW(5) )+8.*GW( 3) )/3«-GMID) 

GAM  ID  =  ABS(GMID)  

GFPS  *  GAMID*GRE+GAE 

IF(GDFLT-GEPS)  1026*  1027  ♦  1027  "" .  ~ 

10  26  IGSWR  =  1 

1027  IF( IGSWA)1029,  1029,  1028 ~ ' 

1028  TF(  IGSWB)  1.029,  1029,  1038 

102^  GHSUG  =  (GEPS/(GFPS+GDELT)+0,45)*GH 

IF( IGSWB) 1030,  1030,  1038  r 

1030  IF(  IGSWA)  1033*  1033»  1031 - -~ 

1031  IF(GH-GHMMI) 1037t  1037*  1032        j- 

1032  IGSWA  =0  

1033  IF(GHSUG-GHMMI )1035»  1034,  1034 

1034  GH  *  GHSUG  


Ill 


GO    TO    1016 

— 103  5-IFJGH-GHMMI  )  1036»    1037  f-1036- ;~       — 

l0?6    GH    =    GHMMI 
•   GO    TO    mi  6 
10"*7    CONTINUE 
h     1038     IF(GAMID-GMAMI » 1040»     1040 #1039 
10  39    GMAMI     =    GAM  ID 

•     -      —    GXMAX    =    GXMOM  •    : 

GHXMA  =  GH  i 

r~ 1040  GRAL  =  GRAL+GH*GMID  \ 

GXMOM  =  GX 

GFXMO  *  GW(5)  

GSIGM  =  GSIGM+GDELT*GH 


•  IF(  IGSWA11012*  1012»  1044 —  — ~ — 

FND 

MON'SS      EXEO  LINKLOAD  —  

CALL  REFLECT 
MCM$$—  EXEO  REFLECT»MJB 
MCN$$      JOB   ACT$$      JIM  ROWLAND     NE         0093U4041 
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9.5.2  SPECTRUM 

The  SPECTRUM  computer  code  was  written  to  apply  the  necessary  correc- 
tion factors  to  the  raw  data  obtained  from  the  neutron  spectrometer.  After 
all  corrections  have  been  applied,  the  data  is  normalized  so  that  the  magni- 
tude of  the  experimental  point  at  the  most  probable  velocity  is  unity. 
Equation  6  is  then  fitted  to  this  normalized  experimental  data  between  se- 
lected points  using  a  least-squares  technique. 

SPECTRUM  consists  of  two  parts,  DATA  and  SUMSQ.  SUMSQ  calculates  the 
squares  of  the  deviations  between  a  trial  theoretical  curve  and  the  normal- 
ized experimental  points.  DATA  computes  everything  else. 

Data  is  read  by   the  program  in  the  following  order: 

1.  A  card  giving  the  number  of  descriptive  cards 

2.  Several  descriptive  cards 

3.  A  card  giving  the  number  of  data  points 

4.  A  card  giving  the  number  of  the  beginning  point  for  the  least- 
squares  fit 

5.  A  card  giving  the  number  of  the  ending  point  for  the  least-squares 
fit 

6.  Cards  giving  the  arm  positions  for  the  data  (as  determined  from 
the  spectrometer  calibration  ring) 

7.  Cards  giving  the  data  channel  counts 

8.  Cards  giving  the  monitor  channel  counts 

9.  A  card  giving  the  number  of  background  points 
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10.  Cards  giving  the  background  arm  angles 

11.  Cards  giving  the  background  data  channel  counts 

12.  Cards  giving  the  background  monitor  channel  counts,  and 

13.  A  card  giving  the  time  taken  for  some  number  of  monitor  counts, 
that  number  of  monitor  counts,  and  the  resolving  time  of  the 
counter. 

Two  different  forms  of  output  are  obtained,  i.e.,  punched  cards  and 
printed  output.  The  punched  cards  are  in  a  form  to  be  used  with  a  suitable 
plotting  program.   This  output  is  in  the  following  form: 

1.  A  card  giving  the  number  of  descriptive  cards 

2.  Several  descriptive  cards 

3.  A  card  giving  the  number  of  data  points,  and 

4.  Cards  giving  velocity,  normalized  experimental  countrate,  energy, 
and  theoretical  countrate. 

The  printed  output  consists  of  the  following: 

1.  A  listing  of  the  descriptive  cards 

2.  The  value  of  the  first  guess  at  K..  in  equation  6 

3.  The  value  of  the  first  guess  at  v  as  given  in  equation  6 

4.  A  table  containing  many  pertinent  values,  each  of  which  is 
identified  by  an  appropriate  table  heading,  and 

5.  The  least-squares  fitted  constants  for  K-  and  v  above. 
The  important  variables  used  in  DATA  are  explained  in  Table  XVII. 

The  execution  times  for  SPECTRUM  are  highly  variable,  being  between 
5  and  20  minutes  per  data  set.  A  long  execution  time  generally  indicates 
a  poor  initial  guess  at  either  v  or  K. . 
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Table  XVII.  Important  Variables  in  DATA 


Symbol 


Explanation 


ACOUNT 

AINC 

AMAX 
ANTER 

ARM 

B 

BARM 

BCOUNT 

BG 

BGM 

BINC 

CHECK 

CONST 

CORR 

DEAD 

EFF 

ENERGY 


an  array  containing  the  corrected  normalized  experi- 
mental data 

the  present  value  of  the  increment  for  changing  K. 

the  normalizing  value  for  K. 

an  array  containing  the  geometrical  correction  factor 
as  described  in  section  2.5.2 

an  array  containing  the  arm  positions  for  the  data 
points 

conversion  factor  from  degrees  to  radians 

an  array  containing  the  background  arm  positions 

an  array  containing  the  theoretical  values  for  the 
leakage  spectrum 

an  array  containing  the  data  channel  counts  for  the 
background 

an  array  containing  the  monitor  channel  counts  for 
the  background 

the  present  value  of  the  increment  in  the  v 

an  array  containing  the  monitor  channel  counts  for 
the  data 

1. 

C(v)  as  defined  in  section  2.7 

the  resolving  time  of  the  data  channel  counting  system 
in  microseconds 

an  array  containing  the  detector  efficiency  for  the 
data  channel  BF~  probe 

neutron  energy 
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Table  XVII.  (continued) 


Symbol 


Explanation 


11 

L 

LA 

MON 
N 

NB 

NR 

REF 

SAINC 

SBINC 

SUM 

SUMA 
THEORY 

THETA 

TIME 

VEL 

WAVE 

XCOUNT 
ZERO 


a  counter  for  determining  the  number  of  iterations  in 
the  curve-fitting  procedure 

the  first  data  point  to  be  fitted  to  the  theoretical 
curve 

the  last  data  point  to  be  fitted  to  the  theoretical 
curve 

the  monitor  counts  for  the  resolving  time  correction 

the  number  of  data  points 

the  number  of  background  points 

the  number  of  descriptive  cards 

the  crystal  reflectivity  as  described  in  equation  16a 

the  initial  value  of  the  K.  increment 

the  initial  value  of  the  v  increment 

the  sum  of  the  square  of  the  deviations  between  the 
theoretical  and  experimental  points 

the  previous  value  for  SUM 

an  array  containing  the  final  ;  eoretical  points  for 
the  spectrum 

an  array  containing  the  arm  angles  for  the  data 

the  resolving  time  of  the  data  channel  counting  system 

an  array  containing  the  neutron  velocities  for  the 
corresponding  arm  angles 

an  array  containing  the  neutron  wavelengths  for  the 
corresponding  arm  angles 

an  array  containing  the  data  channel  counts  for  the  data 

the  zero  point  of  the  spectrometer  calibration  ring 
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M0N1*      JOB   SPECTRUM 
MCN^a,      CGMT  SAVE  PUNCHED  CARDS 
MCNii      CCMT  15  MINUTES»07  PACES 
MONSJ.      ASGN  MJB«12 
MONS*      MODE  GO»TEST 
MCNSS      ASGN  MGC,16 
•  MCW7  1-      EXEQ  FORTRAN,, ♦,,, ,DATA 

RFAL^CN  ' .  '  ""  ... 

C      N  SHOULD  NOT  EXCEED  75 
C      THIS  PROGRAM  GOOD  EOR  200  PLANE  ONLY 

DIMENSlONARM( 75) *VEL( 75) t WAVE  I  75) ♦ ENERGY! 75) »RLFt  75) »E 
.   1FF<75)»ANTER<  ?5 ) »CORR ( 75 )» THEORY ( 75 ) , BCOUNT ( 75 ) » THETAl 75 ) » 
2AC0UNK75  )  ,XCGUNT(75  )  » CHECK*  75) 

DIMENSIONBARM( 75) »BG( 75) »BGM( 75) 

DIMENS10NCHI  (21  ) 

D I  MENS  I ONSUMSOA ( 2  )  ~~ 

CCMMCNVEL* BCOUNT »ACCUNT 

EOU I  VALENCE ( BCOUNT ( 1 ) ,THETA( 1  )  ) 

EOU I  VALENCE ( CHECK (  1  )  »ANTER ( 1  )  ) 

EQUIVALENCE! BARM< 1 ) »VELI1 )>»<  BG( 1)  ,*AVE( 1 )  ) » (BGM( 1 )  ,ENE 
1RGY( 1 ) ) 

1  FORMAT ( 13 ) 

2  FORMAT (5F6.0) 

'3    FORMAT! 1H1,1X,9HARM  ANGLE , 2X , 8H VELOC I f Y » 2X , lOHwAVELENGTH 
1»2X,6HLNERGY,2X,12HREFLECI I  VI 1 Y»3X»VHDt I  .EFF.,2X,12HINIE 
lRCEPTICNt2X»luHCCRRECT ION, 2X ♦ 12HEXPER I  MENTAL ♦ 3 
1X,11HTHE0RET ICAL) 
4  FORMAT  (  3X,7HDEGREES»5X,5HM/SEC»4x»9HANGi>l  RCMo ,4a , 2HEv ,4a  . 
1  ,?HEV»4X,13H=1  AT220uM/S»3X,7HPERCENT,3X,12Hl/2  IN. 
1BEAM»4X»6HFACTCR»4X»11HCCUNTS/M  IN. »4X , 1 1HCCUNTS/MI N. / ) 

5  FORMAT(F9.2,F12«QtF.9..3,FJl1.4»F10.3»F13.3,2F14.5,F12.5»F1 

14,6) 

6  F0RMAT(5F11.4) 

7  FORMAT!  5F9.U  )     ..    _  -  ...   .  ... 

8  FORMAT! 1  HI ,9HCCNST ANT* » E14. 7/ IX , 1 7HMAX I MuM  VELCCI  I  Y*»F6.U) 

9  FORMAT! 1X»79H 

1  -..-,..-._....        .  .- 

10  FORMAT  ( 1H1 ) 

11  FORMAT! 15 ) 

12  F0RMAT(E14.7). " 

13  F0RMAT(4E14.8) 

J  =  (. - 

KA  =  u 
C      M  IS  THE  NUMBER  CF  OBSERVATIONS 
C      LIS  THE  NUMBER  CF  IHE  FIRbl  PC1NF  IC  Fit 
C      LA  IS  THE  NUMBER  CF  THE  LAST  PCINI  IC  FII 

READ(1»1)NR  __     

" WRITE (2,1 )NR  " ~~~ "* 

WRI TE( 3» lu) 

DC22I=1,NR      . 

READ! 1,9) 
WRITE! 3*9 ) 

22  WRITF(2,9)      ■     " - - 

RFAD! 1,1 )N»L,LA 

WRI"TE(2,11)N  T" 

READ! 1,2)ZERC 


READ! 1,2)  ( ARM(  I  )  ,1  =  1  »N) 
READ! 1,2) (XCCUNT( I ) , I=1,N) 
READ! 1,2) (CHECK! I ) ,I=1,N) 

READ(1,1)NB  

READ(1»2) <BARM( 1 ) >I=1»NB) 


. —  — .  —  ■  — 


,r 


READ(1,2) (BG( I  )*I  =  1 ,NB) 
READ(1»2) (BGM( I ) »I=1,NB) 
READ(1»7)TIME»MCN»DEAD- • 
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DTIME=MON*DEAD/TIME    ' 

DC37I=l,NO 

37  B.G(  I  )=BG(  I  )/BGM<  I)        

I=NB   " 

BARM( 1+1) =BARM< I )+BARM( I )-BARM( 1-1 ) 

w9(  ^rU:?6!  n  +  (BG<  n-bG(i-i)  )*<bari*h  i+n -barmi i  m/i barm ( 

1 1  ) -bARM (  I -1 )  ) 

D043I=1*N _       . 

XCCLNTII) =XCCUNT( I )/ CHECK ( I) 
4  3    XCCUNK  I  )=XCCUNT(  1  )*<  l.+XCCUNK  1  )*L>TIME) 
.    DC54I=1»N 
J  =  l 

46  IF(  BARM(  J)-ARM(  I  )  )  51  »47»49  •- 

47  XCCUNTl I )=XCCUNT( I )-BG( J) 
G0T054 

49  J=J  +  1  v  •:•'-  -~ - 

GCTC46 
51  1F( J.EQ.l )GCTC47 

J=J-1 


5  3  XCCUNT(  1  )=XCCUNT(  1  )  - < BG (J+l ) -BG ( J )  )*<ARM(  I  )-BARMU)  )  /  (  BA 

1RM(  J+1)-BARM(  J)  )-BG(  J)    — . :.— 

54  CONTINUE 

AN=(6.02E23)*57.6/( 76, 0*22 .4E03 ) 

6=3.1415926/ 180,  -    -         -  --- 

DC70I=1 ,N 
„ THETA(  I  )  =  ( ARM(  I  )-ZER0)/4. 

SAB  =  5IN(  (  THETA(  I  )  )*B)  "   "  "   ~    ~ 

WAVE( I )=4.ol*SAB 

VFL( I )=982./SAR 

ENERGY( I )=5.04E-G3/(SAB*SAB) 

A=(8.+l3.*EXP(-2.6*.0081/( ENERGY (  I  )**4)  )  )/8. 

REH  (  I  )  =  (  . 5u2 68+3.308 3 E-06*(THET A <  I  )**2*816)  )/VEL( I ) 

EFF(  I  )=l.-EXP(-AN*(3840.E-24)*5U.ti*5QRt  t.u2  5/ENERGM I )  )  ) 
.  AMTER( I )=1.  _ 

IF(SAB.GE..2  5)GCfC69" 

ANTFR( I )=4.*SAB 

6  9  CCRR( I )=REF( 1 )*EFF( I )*CCS(THETA< I )*B)*ANTER( I )*VEL( I )*VE 

1L ( I )*A 
70  ACCUNK I )=XCCUNT( I )/CCRR( I ) 

AMAX  =  ACCUNT(L)  -   " 

DC76J=L»LA 

IF(AMAX-ACCUNT( J) ) 74.76*76 L   _     _ 

74  AKAX=ACOUNT( J) 

K  =  J 
76  CONTINUE  -•■■     -  -  •■   - 

D078I=1 »N 
78  ACOUNT(  I  )=ACOUNTU  )/AMAX 

7  9  CONST =ACOUNT(K)/( ( VEL ( K ) **3 )*EXP (-( 1.225 ) **2 ) ) 

WRI TE(3»12)CCNST*VEL(K) 

.   IK  =  0  .  _ 

SBINC-2. 
CONSTA=CONST 
CONSTB=VEL(K) 
85  AINC=CONSTA/100. 
IZ= IZ+1 

IF( IZ.GE.6)GCTC114 
SAINC=AINC 
.  .  BINC=0.  .  ..„ 
GOT093 


91  BINC  =  SBINC    "  

AINC=0. 

93  JQ  =  o 

94  CALLSU.MSO(CCN'&TA.CaNSTB"fL7C?r»30MT" 

95  CCNSTA=CCNSTA+AINC  ' 
CCNSTB=CCNSTB+BINC-~ — 

_SUMA=SUM 


116 


CALLSLM5U(CGNbTA»CGNSIBtLtLAtSUM) 
1  I-  (  SUMA.E0.SUM)GGTGlU6 
I F( SUMA.GI .SUM) GOT 09 5 
l£l  I F I JO.EO.l JGCTC106      

BINC=-BINC 

J0=1 

G  0  T  0  9  5 
106  I F( AINC.EQ.O. )G0TOlO8 

G0T091 

1F( AB5(AINC)-(SAINC/10C.) )37t37t48 
108  I  F  (  dINC»E(J.u.  )GGTGllU 

"  GOT 08 5  - -  - 

1FI ABS(BINC)-(SBINC/100. ) )36t36t48 
11U  AINC=-AiNC/lU.  . . 

BINC=-BINC/10. 

J0=0 

G0TG95 

114  1=1 

115  IF( I .J0.1 )GCTC123 

IF(  I  .E0.72  JGCT0123 " •"  ~  " 

117  THEORY ( I )=CGNSTA*(VEL( I ) **3 ) *EXP ( - ( 1 .22  5*VEL ( I J/CCNSTB)* 
1*2)  ) 

IF(  I  •GT.ISDGCTC126 

*RlIE(2tl3)VLL(I) »ACGUNT( I ) *LNEkGY( I ) ♦ THEORY { I ) 

WR1  rE(3»5 )ARM(  1  ) »VtL( 1  )  tWAVEt 1 ) tENLRGY< I ) tRtF<  1 ) tEFF ( I ) ♦ 
lANTERt I )  »CCRR(  I  ) »ACCUNT( I )  ♦  THEORY (  I ) 

1  =  1+1 

GCTC115 
123  WRITE(3*3) 

WRI  TF(3*4)  -- 

GCT0117 
126  CONTINUE 

WRI TF(3»8)C0NSTAtC0NSTB 

CALLEXIT 

END 
MON-t-i      EXEQ  FORTRAN* ttttVtSUMSCT         " 

SUdkOUTlNESUMSQ(CGNSTAtCGNSTBtLtLAtSUM) 

COMMONVEL ( 75 ) tBCOUNT ( 75 ) t ACOUNT (75 ) 
1  F0RMATOE14.7) 

SUM=0# 

D05I=LtLA 

BCOUNT( I ) =CONSTA*(VEL( I )**3 >*EXP ( - ( 1 • 22 5* (VEL( I ) ) /CONSTB 

$  suM=sufo+(6couNT(  i  i-AcouNrrni^i/ACcuNnn 

RFTURN 

END   -  -  -  - 

MGNiS      EXEQ  LINKLOAD 

CALL  DATA 
MONSi'     EXEQ  DATAtMJB   ~ —    


...  -.._ 
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9.5.3  CROSS  SECTION 

The  CROSS  SECTION  computer  code  was  written  to  process  data  obtained 
from  beryllium  total  neutron  cross  section  measurements.  This  program  is 
written  for  measurements  taken  using  beryllium  discs  and  a  Li F  (111)  crystal. 
In  order  to  determine  cross  sections  of  other  materials  or  to  use  other 
diffraction  crystals,  some  of  the  internal  constants  in  the  program  must 
be  changed. 

The  program  is  rather  straight  forward.  Down  to,  and  including,  state- 
ment 38,  the  data  is  merely  read  into  the  computer  and  the  background  sub- 
tracted. 

The  data  input  is  similar  in  nature  to  that  in  the  SPECTRUM  computer 
code.  The  input  is  given  in  the  following  form: 

1.  A  card  giving  the  total  number  of  data  points 

2.  A  card  giving  the  location  of  the  zero  point  as  determined  from 
the  spectrometer  calibration  ring 

3.  Cards  giving  the  arm  positions  for  the  data  as  indicated  by  the 
spectrometer  calibration  ring 

4.  Cards  giving  the  data  channel  counts  for  the  beryllium  attenuated 
neutron  beam 

5.  Cards  giving  the  monitor  channel  counts  for  the  beryllium  atten- 
uated neutron  beam 

6.  Cards  giving  the  data  channel  counts  for  the  unattenuated  beam 

7.  Cards  giving  the  monitor  channel  counts  for  the  unattenuated  beam 
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8.  Cards  giving  the  thickness  of  the  beryllium  discs  (in  numbers 
signifying  the  total  thickness  in  eighths  of  an  inch)  for  each 
different  arm  position 

9.  A  card  giving  the  number  of  background  data  points 

10.  Cards  giving  the  arm  positions  for  the  background  as  indicated  on 
the  spectrometer  calibration  ring 

11.  Cards  giving  the  data  channel  counts  for  the  background,  and 

12.  Cards  giving  the  monitor  channel  counts  for  the  background. 

The  output  consists  of  punched  cards  which  are  suitable  for  use  with 
PLOT  (see  section  8.5.2)  and  a  printed  listing.  The  punched  cards  come  out 
in  the  following  form: 

1.  A  card  giving  the  number  of  data  points,  and 

2.  Cards  giving  the  computed  value  for  the  total  neutron  cross  section 
and  the  corresponding  velocity. 

The  printed  output  consists  of  a  table  (with  column  headings)  listing 
arm  position,  corresponding  neutron  velocity,  wavelength,  and  energy,  and 
the  total  neutron  cross  section. 

The  important  variables  contained  within  CROSS  SECTION  are  summarized 
in  Table  XVIII. 
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'able. XVIII.  Important  Variables  Contained  in  CROSS  SECTION 


Symbol 


Explanation 


ACOUNT 

ARM 

AZERO 

BARM 

BG 

BGM 

DATA 

DMON 

ENER 

N 

NB 

SIGMA 

VEL 

WAVE 

XCOUNT 


an  array  containing  normalized  unattenuated  beam  data 

an  array  containing  arm  positions  for  data 

the  zero  position  of  arm 

an  array  containing  arm  positions  for  background 

an  array  containing  data  channel  counts  for  background 

an  array  containing  monitor  channel  counts  for  back- 
ground 

an  array  containing  data  channel  counts  for  Be 
attenuated  beam 

an  array  containing  monitor  channel  counts  for  Be 
attenuated  beam 

neutron  energy  in  eV 

number  of  data  points 

the  number  of  background  points 

the  total  neutron  cross  section  in  barns 

neutron  velocity  in  meters  per  second 

neutron  wavelength  in  angstroms 

an  array  containing  normalized  attenuated  beam  data 
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M  ON' LI- 

JOB 

MCNSS,, 

CCMT 

MO\'$$ 

ASGN 

MONT-S 

ASGN 

MCMXJ 

MODE 

MONiS 

EXEU 

CROSS  SECTION 

15»5»ROWLAND,NE 

MGC , ] 6 

MJR,12 

c  ^ 

FORI  RAN,  ,  »»*  ,  ,R0ULANL> 
. ..   OIMENSlONUATA(  10U  )  ,OMON(  luO)»RDATA(loO)  ,ROMON(  100  ) 
l»THICK(luu),ARM(luo)  »XCOUN  F  (  1  00  )  »  d<j  <  1  uo  )  ♦  bGM  (  luu  ) 
0IMENSI0NBARM(  100)  »BG(  100)  »BGM(  luu  ) 

4  FORMAT (  I  3 ) 

5  FORMAT ( 5F6.u) 

6  FORMAT! 1  HI ,  1 X  ,  9HARM ANGLE » 2X , 8HVEL0C I TY ♦ 2X , 1 OHw AVE  LENGTH 
l»2X»6HENERGY»5X»5HSIGMA) 

8  FORMAT ( 5X, 3HCM •♦ 7X, 5HM/ SEC »4X»9H ANGSTROMS ,4Xf2HEVf8Xf 
..15HRARNS) 
10  FORMAT(F9.2,Fl2.0,F9.3,FlTV47frlO.3j  " 

1  F0RMAT(2E14.8) 

TOTAL  CROSS  SECTION  FOR  BERYLLIUM  USING  LIF  (111)  CRYSTAL 

READ(1,4)N 

READ( 1 ,5 ) AZERO 

REAO( 1,5) (ARM( I ) ,1=1 »N)   - 

READ( 1 ,5) (DATA( I ) , 1=1 »N) 

REA0( i,5) (DMONI I  )  , 1  =  1 ,N) 

RFAD( 1,5)  (RDATA(  I  )  »  1  =  1  ,Nf 

RFAD< lf 5) (RDMON( I ) » 1=1 ,N) 

RFAD( 1,5)  (THICK(  I ) »  1=1 »N) 

RFAD( 1 ,4) NB 

RFAD( 1 ,5) <BARM( I ) , 1=1 ,NB) 

READ( 1,5) (BG( I ) ♦ 1=1 tNB)        

READ( lfb) (BGM( I ) »I=1»NB) 

D016I  =  1,N  

"  '  XCOUNTI  I  )=DATA{  n/DMONTT) 

16  ACOUNT( I ) =RDATA( I )/RDMON( 1 ) 

0018  1  =  1  fNB  -  — * ... 

18  BG(  I  )=BG(  I  )/BGM(  I ) 

C0NST=8./( . 1236*2.54) 

B=3. 1415926/ 180, 

WRITE(3f6) 

WRITE!  3,8)  .  .      _     __ 

.  _  WRI TE(2,4)N 

BARM(NB  +  1)=BARM(NB)*2.-BARM<N'B-1)  '•" 

bG(NB+l )=BG(NB)+(BG(NB)-BG(NB-1 ) )*(BARM(NB+1 )-BARM(NS) )/ 
1 (BARM(NB)-BARM(NB-l) ) 
"   00381=1, N   '      

J  =  l 

28  IF(bARM( J)-ARM( I )) 34,29,32   -  -   -■  

29  XCOUNK I )=XCOUNT( I )-BG( J) 
ACCUNT< I )=ACOUNT( I )-BG( J) 
G0T038  - - 

32  J=J+1 
G0T028 

34  IF( J.EQ.l JG0T029  "" 

J  =  J-1 
3  6  XCOLNT( I  )=XCOUNT( I  ) - ( BG( J+l )-BG( J )  )*(ARM(  I )-BARM( J)  )/ 
1(BARM( J+l )-BARM( J) )-BG( J) 
ACOONTt I )=ACOUNT( I ) - ( BG ( J+ 1 ) -BG ( J ) )*(ARM( I )-BARM( J) )/( 

--•1BARM(  J+1)-BARM(J)  )-BG<  J)  - -  -  r-    - 

38    CONTINUE 

D046I=1,N  •' 

~~  SIGMA=(CONST/THICK(  D")  *ALCG(  ACCUNTtT)  AXCCUNTTI  )  7 r^" 

SAR=SIN( ( ARM( I ) -AZERO )*B/4. ) 

..    WAVE  =  4.64*SA8     ......  -..- -—.-    -  -   -      - '-■ 

VFL=851./SAB 
EMER=3.78E-0  3/(SAB*SAB) 
--  WRT  TF(  7  «  }  )  VFl-fSTGMA      ■         ■ — '--- :-      -   ■  ■  —  •.-.-—  .-  — 


-  i 
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46  WRI  TE(3»lu)ARM( I ) » VEL ♦ WAVE »ENER »S IGMA 
....._  CALLER I  T : 

MCNiS  .     EXEO  LINKLCAD 

.--- CALL  ROWLAND    .•—...,..... 

MCNSS      EXEQ  RCWLANDjMJB 
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.     Determination  of  an  Empirical  Equation 
for  Integrated  Reflectivity 

The  solutions  to  equation  15  .  :ed  by  the  R  FIVITY  computer 

program  (see  Appendix  E  for  explanation  of  this  program)  are  given  in  Table 
XIX. 

As  indicated  in  section  2.5.1,  let  us  try  to  fit  the  above  data  in  Table 
XIX  with  an  empirical  equation  of  the  following  form: 

ft 

—    =  D  +  K3k.  (16a) 

sine 

In  order  to  eva'.uate  the  three  empirical  constants,  D,  H,  and  k,  thi 
different  values  of  6  were  chosen,  i.e.,  15°,  35°,  and  50°.  Thus,  three 
simultaneous  equations  in  three  unknowns  were  obtained  which  were  then 
solved  by   conventional  methods.  The  values  for  the  constants  found  in  this 
>r  are  given  in  Table  XX. 

The  approximation  is  thus  accurate  to  about  3%   of  the  desired  va". ue 
for  both  planes  and  for  the  range  of  the  tables,  and  accurate  to  about  1/2% 
within  the  range  of  the  data  taken  in  these  experiments. 
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table  XIX.  Summary  of  Results  from  REFLECTIVITY  Computer  Progr 


Diffraction 

A 

R  (minutes 

■  • J 

)  A 

Re/sine 

(minutes  )x  10 

angle 

111  plane 

200  plane 

111  plane 

200  plane 

5 

.10304 

.043841 

1.1822 

.50299 

10 

.20501 

.087542 

1.1805 

.50413 

15 

.30754 

.13179 

1.1882 

.50920 

20 

.41231 

.17732 

1.2055 

.51845 

25 

.52115 

.22494 

1.2331 

.53225 

30 

.63618 

.27559 

1.2724 

.55118 

35 

.75997 

.33047 

1.3250 

.57615 

40 

.89585 

.39113 

1.3937 

.60849 

45 

1.0482 

.45967 

1.4824 

.65007 

50 

1.2235 

.53910 

1.5972 

.70375 

55 

1.4309 

.63390 

1.7468 

.77385 

60 

1.6849 

.75114 

1.9456 

.86736 

Table  XX.  Summary  of  Empirical  Constants  for  Equation  16a 


Parameter 


Crystal   plane 


: : " 


Oi^ft 


1.1755  x  10'3  min.     .50268  x  10"3  min, 


1.9898  x  10~9  min.  3.3033  x  10"9 


2.907 


2.C.6 


L       .   .   I 
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T  ACT 

A  neutron  <  :ed  collimator  have  been  con- 

structed by  the  author.  This  equipment  has    1  used  to  study  the  neutron 

a  from  the  d'  the  KSUTMII  reactor, 

effects  upon  these  spectra  of  changing  tl    .1  loading  pattern 
been  s     .It     jwn  that  for  certain  types  of  experiments,  proper 
selection  of  the  beamport  and  fuel  loading  is  important  in  obtaining  the 
)0ssible  signal  to  noise  ratio. 

s  rements  have  also  bee.,  made  of  the  total  neutron  cross  s_ 
for  metallic  beryllium  between  the  energies  of  .0057  and  .045  eV.  Higher 

order  diffractions  c the  results  to  be  in  error  by   as  much  as  300%  below 

.  lz   .CC53  eV.  Above  ....  eVs  the  results  are  generally  within  2. 
the  published  values.  The  break  in  the  cross  £^c;;c:.  at  .021  eV  "is  yary 
.'?,  indicating  a  resolution  of  .006  >_ .  at  :his  e 

:.as  been  de    -crated  that  this  neutron  spectrometer  may  be  oper- 
ted  safely  provided  that  the  instructions  in  section  3.4  are  closely  fol- 

: In  this  connection,  the  author's  total  maximum  film  badge  reading 

j   calendar  quarter  was  ISO  mrem. 


